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A STRONG FORM OF THE QUANTITATIVE WULFF INEQUALITY 


ROBIN NEUMAYER 


Abstract. Quantitative isoperimetric inequalities are shown for anisotropic surface energies where the 
isoperimetric deficit controls both the Fraenkel asymmetry and a measure of the oscillation of the boundary 
with respect to the boundary of the corresponding Wulff shape. 


1. Introduction 

1.1. The Wulff inequality and stability. For n > 2, the well-known isoperimetric inequality states that 

P{E) > n\B\ 1/n \E\ 1/n ', 

n' = n/(n — 1), with equality if and only if E is a translation or dilation of B = {x £ R ra : |x| < 1}, the 
Euclidean unit ball in R n . The isoperimetric inequality holds for all sets of finite perimeter E C I", with 
the perimeter of E equal to P{E) = 'H n ~ l {d*E). Here, d*E is the reduced boundary of E; see Section 
or f3(Tj for a more complete overview. 

The anisotropic surface energy is a natural generalization of the notion of perimeter and has applications 
in modeling of equilibrium configurations for solid crystals (see [33] [27], [33]) and of phase transitions (see 
[23j). We introduce a surface tension f : R n —> [0, +oo) to be a convex positively 1-homogeneous function 
that is positive on 5 ra_1 . The corresponding (anisotropic) surface energy of a set of finite perimeter E C R" 
is defined by 

$(£)= [ f{v E {x))dU n -\x). 

J d* E 

Here, v E is the measure theoretic outer unit normal; see Section [2. 1[ Just as the ball minimizes perimeter 
among sets of the same volume, as expressed by the isoperimetric inequality, the surface energy is uniquely 
minimized among sets of a given volume by translations and dilations of a fixed convex set K determined 
by the surface tension. This set K is given by 

K = Q {x £ R" : x ■ v < /(i/)}, 

veS"- 1 

and is known as the Wulff shape of $. The minimality of the Wulff shape is expressed by the Wulff inequality: 

*(E) > W(j^) 7 = n\K\V n \E\ l / n ’, 

with equality if and only if E is a translation or dilation of K ; see [Mi mi ua [sni eh uni ni- In the case 
where the surface tension / is constantly equal to one, $(E) reduces to the perimeter P(E) and the Wulff 
inequality reduces to the isoperimetric inequality. 

Given a surface tension /, one may define the gauge function /* : R" —> [0, +oo) by 

f*{x) = sup{x • v : f{v) < 1}. 

The gauge function provides another characterization of the Wulff shape: K = {x : f*(x) < 1}. 


2.1 


l 
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To quantify how far a set is from achieving equality in the Wulff inequality, we introduce the anisotropic 
isoperimetric deficit , or simply the deficit , of a set E , defined by 


<M#) = 


3>(U) 




- 1. 


The deficit equals zero if and only if, up to a set of measure zero, E = x + rl\ for some x € R n and r > 0. 
This quantity is invariant under translations and dilations of E, as well as modifications of E by sets of 
measure zero. 

The optimality of the Wulff shape in the Wulff inequality naturally gives rise to the question of stability: 
does the deficit control the distance of a set from the Wulff shape? In other words, given a distance d from 
the family {x + rK : x £ r > 0}, one wants to find inequalities of the form 


6*(E)>w(d(E)), 


( 1 . 1 ) 


where w is a (possibly explicit) function such that oj(d{E)) -A 0 + as d(E) —> 0 + . Ideally, one hopes to find 
the function uj that provides the sharp rate of decay. Such an inequality can be viewed as a quantitative 


form of the Wulff inequality: by rearranging the deficit, (1.1) becomes 


HE) > *(K) ( j|| 


i A 


+ w{d{E))n\K\ 1 / Tl \E\ 1 / n '. 


In this way, iv(d(E)) serves as a remainder term in the Wulff inequality. 
A well studied distance is the asymmetry index , a$(E), defined by 


. ( \EA(rK + y)\ 

<**&)= mm |- — - 


\rK\ = \E\ 


( 1 . 2 ) 


where EAF = (E \ F) U (F \ E) is the symmetric difference of E and F. For the case / constantly 
equal to one, the asymmetry index is known as the Fraenkel asymmetry. The quantitative isoperimetric 
inequality with the Fraenkel asymmetry was proven in sharp form by Fusco, Maggi, and Pratelli in [23]. 
Using symmetrization techniques, they showed that if If is a set of finite perimeter with 0 < \E\ < oo, then 


ai(E) 2 < C(ra)<5i(£7). 


(1.3) 


Here and in the future, we use the notation <5i and oq for the deficit and asymmetry index corresponding to 
the perimeter. 


Before this full proof of (1.31 was given, several partial results were shown; see ED [23 [2g . Another proof 
of (1.3) was given in [8], introducing a technique known as the selection principle, where a penalization 
technique and the regularity theory for almost-minimizers of perimeter reduce the problem to the case 
shown in [21] . 

Stability of the Wulff inequality was first addressed in m, without the sharp exponent. Figalli, Maggi, 
and Pratelli later proved the sharp analogue of (1.3) for the Wulff inequality in [IS], using techniques from 


optimal transport and Gromov’s proof of the Wulff inequality. They showed that there exists a constant 
C(n), independent of /, such that 

a$(U) 2 < C(n)6*(E) (1.4) 


for any set of finite perimeter E with 0 < \E\ < oo. In both (1.3) and (1.4), the power 2 is sharp. The 
constant C(n) in (1.4) is explicit, a feature that is not shared with any other proofs of (1.3), (1.5), or any 
of the quantitative inequalities obtained in this paper. 

In [22], Fusco and Julin proved a strong form of the quantitative isoperimetric inequality, improving (1.3) 
by showing 


ai [Ef+^Ef < C{n)Si{E) 


(1.5) 
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for any set of finite perimeter E with 0 < \E\ < oo, where the oscillation index (3i(E) is defined by 


Pi(E) = min 


yeR" \ \2n\B\ 1 / n \E\ 1 / n ' J d , E 


r 2 \ 1/2 ■n 

J ve{x) - VB r {y) (y + _ y 'j dU n ~ 1 (x)\ : \B r \ = \E\ j 

1/2 


= mm 

yeR n V n\B 


V n \E I 1 /”' Le 


1 t~ - ■v E (x)d'H n 1 (a;) 

\ x - y I 


( 1 . 6 ) 


While the asymmetry index ot\(E) is an L 1 distance between E and B , the oscillation index /3\(E) quantifies 
the oscillation of d*E with respect to dB. In [25J Proposition 1.2], Pi(E) is shown to control a\(E)\ see 
Proposition |2.4| for the analogous statement in the anisotropic case. Once again, the power 2 in is 
sharp for both ot\{E) and 

Since [22] . analogous strong form quantitative inequalities have been studied in several settings: in Gauss 
space Piili 3], on the sphere [4], and in hyperbolic n-space [5]. 


1.2. Statements of the main theorems. The goal of this paper is to address the question of what form 
the inequality (1.5) takes in the case of the anisotropic surface tension $. In other words, we prove a 
strong form of the quantitative Wulff inequality, improving (1.4) by adding a term to the left hand side that 
quantifies the oscillation of d*E with respect to dK. We define the <b -oscillation index by 


@&(E) = min 

y6R™ 


1 


/ /{ve{x)) - v E {x) ■ — V —d'H nl {x) 

I d*E f*[x-y ) 


1/2 


(1.7) 


The following theorem is a strong form of the quantitative Wulff inequality that holds for an arbitrary surface 
energy. 

Theorem 1.1. There exists a constant C depending only on n such that 

at,(E) 2 + /3<i,(£0 4n/(n+1) < C6<s,(E) (1.8) 

for every set of finite perimeter E with 0 < \E\ < oo. 


As in (1.4), the constant is independent of /. We expect that, as in (1.5), the sharp exponent for j3$(E) 
in (1.8) should be 2. With additional assumptions on the surface tension /, we prove the stability inequality 
in sharp form for two special cases. 

Definition 1.2. A surface tension f is X-elliptic, A > 0, if f £ C 2 ( 


(\7 2 f(v)r)- t> 


r — 


( t -r) 


\ {0}) and 
2 


for u, r £ K" with v ^ 0. 

This is a uniform ellipticity assumption for V 2 /(u) in the tangential directions to v. If / is A-elliptic, 
then the corresponding Wulff shape K is of class C 2 and uniformly convex (see [22], page 111). When is 
a surface energy corresponding to a A-elliptic surface tension, the following sharp result holds. The constant 
depends on w$ and M$, a pair of constants defined in (2.2) that describe how much / stretches and shrinks 
unit-length vectors. 


Theorem 1.3. Suppose f is a X-elliptic surface tension with corresponding surface energy <E>. There exists 
a constant C depending on n, A,m<j>/M$, and || "X7 2 f\\c°{dK) such that 

a^(E) 2 + /3$(E) 2 < C6$(E) (1.9) 

for any set of finite perimeter E with 0 < \E\ < oo. 

The second case where we obtain the strong form quantitative Wulff inequality with the sharp power is 
the case of a crystalline surface tension. 
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Definition 1.4. A surface tension f is crystalline if it is the maximum of finitely many linear functions, 
in other words, if there exists a finite set {xj}J =l C R" \ {0}, N £ N, such that 

f(v) = max {xj • v} for all v e S'"” 1 . 
i <j<N l 0 J J 

If / is a crystalline surface tension, then the corresponding WulfF shape K is a convex polyhedron. In 
dimension two, when / is a crystalline surface tension, we prove the following sharp quantitative Wulff 
inequality. 

Theorem 1.5. Let n = 2 and suppose f is a crystalline surface tension with corresponding surface energy 
$. There exists a constant C depending on f such that 

a$>(E) 2 + /3®(E) 2 < C5^{E) 
for any set of finite perimeter E with 0 < \E\ < oo. 

Some remarks about the definition of the <E>-oscillation index /3$ in are in order. The oscillation 

index (3i(E) in (1.6) measures oscillation of the reduced boundary of a set E with respect to the boundary 
of the ball. Indeed, the quantity 0\(E ) is the integral over d*E of the Cauchy-Schwarz deficit 1 — ^ • v E {x), 
which quantifies in a Euclidean sense how closely veM aligns with At. 

To understand ( |1.7| ), we remark that / and /* are dual in the sense that they yield a Cauchy-Schwarz-type 
inequality called the Fenchel inequality, which states that 

X 

”e(x) ■ YJx) - 

Just as the oscillation index f3\(E) quantifies the overall Cauchy-Schwarz deficit between ^ and ve{x), the 


term /3$(E) is an integral along d*E of the deficit in the Fenchel inequality. In Section 2.2 we show that 
f(yE{x)) = y ■ ve(x) for y £ dK if and only if y is a point on dK where ve{x) is normal to a supporting 
hyperplane of K at y. In this way, (3®(E ) quantifies how much normal vectors of E align with corresponding 
normal vectors of K , and therefore provides a measure of the oscillation of the reduced boundary of E with 
respect to the boundary of K. Note that in the case / constantly equal to one, /?$ agrees with /3i. 


It is not immediately clear that (1.7) is the appropriate analogue of (1.6) in the anisotropic case. Noting 
that x i—^ {x — y)/f*{x — y) is the radial projection of R" \ {0} onto dK + y, one may initially want to 
consider the term 

1/2 


P%{E) = min 
y£R n 


= mm 
y£R n 


1 


2n\K\ 1 / n \E\ 1 / r 

1 


7 [ v E {x) - v K ( 77 — V -^) d'H n 
Jd*E \f*(x-y)/ J 

[ l-v E {x)-v K (- f X ^—) d'H n ~ l {x)\ 

J d* e \f*( x -y) J J 


( 1 . 10 ) 


K n\K\yn\E\x/n Jd , E 

However, in Section [ 6 ] we see that such a term does not admit any stability result for general $. Indeed, in 
Example 6.1 we construct a sequence of crystalline surface tensions that show that there does not exist a 
power er such that 

PUE) a <C(n,f)6*(E) (1.11) 


for all sets E of finite perimeter with 0 < \E\ < oo and for all <!>. Furthermore, Example 6.2 shows that 
even if we restrict our attention to surface energies which are 7 -A convex, a weaker notion of A-ellipticity 
introduced in Definition 1.6 an inequality of the form (1.11) cannot hold with an exponent less than a = 4. 
The examples in Section 6 ] illustrate the fact that, in the anisotropic case, measuring the alignment of normal 
vectors in a Euclidean sense is not suitable for obtaining a stability inequality for general <I>; it is essential 


to account for the anisotropy in this measurement. The <h-oscillation index /3<p(E) in (1.7) does exactly this. 


In the positive direction, when the surface tension / is 7 -A convex, /3 $(E) is controlled by (3<s>(E ). As one 


expects from Example 6.2, the exponent in this bound depends on the 7 -A convexity of /. We now define 
7 -A convexity. 
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Definition 1.6. Let f : R n —► R be a nonnegative, convex, positively one-homogeneous function. Then we 
say that f is 7 -A convex for 7 > 0, A > 0 if 

2+7 


f{u + t) + /(u - t) - 2 f(v) > 


( T 'Ui) 


( 1 . 12 ) 


for all v, t £ R n such that y / 0. 


Dividing (1.12) by r 2 , the left hand side gives a second difference quotient of /. While A-ellipticity 
assumes that / € C 2 (R rl \ {0}) and that its second derivatives in directions r that are orthogonal to v are 
bounded from below, 7 -A convexity only assumes that the second difference quotients in these directions 
have a bound from below that degenerates as r goes to 0. Of course, a 0-A convex surface tension / with 
/ £ C 2 (R ra \ {0}) is A-elliptic. The I v norms f v (x) = CCr=i \ x i\ P Y^ p f° r P € (l,oo) are examples of 7 -A 
convex surface tensions; see Section [6j When / is a 7 -A convex surface tension, the following theorem shows 
that /3$ controls /3|. 

Theorem 1.7. Let f be a 7 -A convex surface tension. Then there exists a constant C depending on 7 , A, 
and /M$ such that 


P 1(E) 


( 2 + 7)/2 


< c 


P(E) 


for any set of finite perimeter E with 0 < \E\ < 00. 


7/4 


/?$(£)• 


As in Theorem 1.3 the constant depends on and which are defined in (2.2). As an immediate 
consequence of Theorem |1.7| Theorem 0 and Theorem [OJ we have the following result. 

Corollary 1.8. If f is a 7 -A convex surface tension, then there exists a constant C depending on 71 , 7 , A, 
and /M§ such that 

7n/(n+l) 


a*(E) 2 + Pl(EY < C 


P(E) 


8*{E) 


^n\K\ 1 / n \E\ 1 / n ' _ 

for any set of finite perimeter E with 0 < |i?| < oo, where a = 2n(2 + 7 )/(n + 1). 

If f is a \-elliptic surface tension, then there exists a constant C depending on n, 7 , A, and 

||V 2 /|| C c (aif) such that 

a*(E) 2 + Pl(E) 2 <C6*(E) 
for any set of finite perimeter E with 0 < |2£| < 00. 


1.3. Discussion of the proofs. At the core of the proof of (1.5) are a selection principle argument, the 
regularity theory of almost-minimizers of perimeter, and an analysis of the second variation of perimeter. 
Indeed, with a selection principle argument in the spirit of the proof of ( |1.3| ) by Cicalese and Leonardi 
in [Hj, Fusco and Julin reduce to a sequence {Fj} such that each F :] is a (A,ro)-minimizer of perimeter 
(Definition 4.4) and Fj —> B in L 1 . Then, by the standard regularity theory, each set Fj has boundary given 
by a small C 1 perturbation of the boundary of the ball. This case is handled by a theorem of Fuglede in [21], 
which says the following: Let E be a nearly spherical set, i.e., a set with barycenter bar E = |1 ?| _1 f E xdx 
at the origin such that \E\ = \B\ and 

dE = {a: + u(x)x : x £ dB} 

: with u £ C x (dB). There exist C and e depending on n such that if HuHcqaB) < e, then 

IMIndas) — C5\(E). (1-13) 


for u:dB 


The proof of (1.13) makes explicit use of spherical harmonics to provide a lower bound for the second 
variation of perimeter. It is then easily shown that aq (E) + j3\(E) < C\\u\\H^(dB), and therefore (1.13) 
implies (1.5) in the case of nearly spherical sets. Indeed, aq (E) < Cfd\(E) as shown in Proposition |2.4[ 
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and in the case of nearly spherical sets, the oscillation index j3\ is essentially an L 2 distance of gradients: if 
w(x) = x + u(x)x , then 

. x(l + u(x)) + X7u(x) 

ve{w(x )) = — , , 

V(l+«) 2 + |Vu| 2 

where the Vu is the tangential gradient of u. Then 

n\K\/3i(E) 2 < [ l-iy E {w)~dH n ~ 1 = f y/{\ + u ) 2 + |Vu| 2 - (1 + u) dU n ~ l 

JOE M Job 

= [ hvu\ 2 + 0(|Vu| 2 ) dU n ~ x < \\u\\ 2 HHdB) . 

Job 1 


In each of Theorems 00 and |1.5[ at least one of the three key ingredients of the proof of Fusco and 
Julin is missing. The proof of Theorem 0 uses a selection principle to reduce to a sequence of (A ,r 0 )~ 
minimizers of $ converging in L 1 to K. However, for an arbitrary surface tension, uniform density estimates 
(Lem ma |3.3| ) are the strongest regularity property that one can hope to extract. We pair these estimates 
with (1.4) to obtain the result. 

The proof of Theorem 1.3 follows a strategy similar to that of the proof of (1.5) in [22]. If / is a A- 


elliptic surface tension, then (A, ro)-minimizers of the corresponding surface energy 4> enjoy strong regularity 
properties. Using a selection principle argument and the regularity theory, we reduce to the case where dE is 
a small C 1 perturbation of dK. The difficulty arises, however, in showing the following analogue of Fuglede’s 


result (1.13) in the setting of the anisotropic surface energy. 


and m$/M$ such that if 
(1.14) 


Proposition 1.9. Let f be a X-elliptic surface tension with corresponding surface energy $ and Wulff shape 
K. Let E be a set such that \E\ = \K\ and bar FI = bar AT, where bar E = lAj -1 :dx denotes the 
barycenter of E. Suppose 

dE = {ac + u(x)vk{x) : x £ dK} 

where u : dK —► M is in C 1 (dIL). There exist C and ei depending on n, A 
I Ml CHaic) < e i) then 

IMIffqdic) — C6$(E). 

Again, ra$ and are defined in (2.2). To prove ( 1.13| ), Fuglede shows that, due to the volume and 
barycenter constraints respectively, the function u is orthogonal to the first and second eigenspaces of the 
Laplace operator on the sphere. This implies that, thanks to a gap in the spectrum of this operator, functions 
satisfying these constraints satisfy a Poincare inequality with a larger constant than the Poincare inequality 
that holds for w £ H 1 (dB) with mean zero (i.e., satisfying only the volume constraint). Fuglede’s reasoning 
uses that fact that the eigenvalues and eigenfunctions of the Laplacian on the sphere are explicitly known. 

The analogous operator on dK arising in the second variation of $ also has a discrete spectrum, but one 


cannot expect to understand its spectrum explicitly. Instead, to prove (1.9), we exploit (1.4) in order to 
obtain a Poincare inequality with a larger constant for functions u £ H 1 (dK) satisfying the volume and 
barycenter constraints. 

Then, as in the isotropic case, one shows that a$(E) + /?$(!?) < Cj|u[|ffi(3R:) f° r a constant C = 
f\\c°(dK)), and therefore (1.14) implies (1.9) for small C 1 perturbations. Indeed, Proposition 


2.4 


C( _ _ 

implies that a$(E) < C(n)/3q,(E), and the fact that (3$(E) < C\\u\\H^(dK) is a consequence of a Taylor 


expansion and a change of coordinates. The computation is postponed until (4.16) as it relies on notation 
introduced in Section [4j 


The proof of Theorem 1 1.5 1 also uses a selection principle-type argument to reduce to a sequence of almost- 
minimizers of <F converging in L 1 to the Wulff shape. In this case, a rigidity result of Figalli and Maggi in 
m allows us reduce to the case where E is a convex polygon whose set of normal vectors is equal to the 
set of normal vectors of K. From here, an explicit computation (Proposition |5 . lj) shows the result. 
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The paper is organized as follows. In Section [2j we introduce some necessary preliminaries for our main 
objects of study. Section [3] is dedicated to the proof of Theorem o while in Sections [4] and [5| we prove 
Theorems 1.3 and 1.5 respectively. In Section [dj we consider the term f3$(E) defined in ( 1.10| , providing 
two examples that show that one cannot expect stability with a power independent of the regularity of / 


and proving Theorem 1.7 


Acknowledgments: The author would like to thank Alessio Figalli and Francesco Maggi for their mentor¬ 
ship, guidance, and many helpful discussions. Further thanks are due to a thorough referee for providing 
several useful remarks. This research was supported by the NSF Graduate Research Fellowship under Grant 


No. DGE-1110007. 


2. Preliminaries 

Let us introduce a few key properties about sets of finite perimeter, the anisotropic surface energy, and 
the dJ-oscillation index /3$. 

2.1. Sets of finite perimeter. Given an M"-valued Borel measure y on M", the total variation |/i| of p, on 
a Borel set E is defined by 

ME) = sup | X] HEj)\ : Ej n£, = 0, [j E 3 C E 
^jeN jew 

A measurable set E C R" is said to be a set of finite perimeter if the distributional gradient D xe of the 
characteristic function of E is an Revalued Borel measure on R" with |D\£;|(R") < oo. 

For a set of finite perimeter E, the reduced boundary d*E is the set of points x £ I" such that 
\T>XE\(B r (x)) > 0 for all r > 0 and 

lim exists and belongs to S'" -1 . (2.1) 

r->0+ \DxE\(B r (x)) 

If x £ d*E , then we let —v E denote the limit in <& We then call v E : d*E —> S n 1 the measure theoretic 
outer unit normal to E. Up to modifying E on a set of Lebesgue measure zero, one may assume that the 
topological boundary dE is the closure of the reduced boundary d*E. For the remainder of the paper, we 
make this assumption. 


2.2. The surface tension and the gauge function. Throughout the paper, we let 


= inf /(u), 

i/ES 1 "- 1 


= sup f(y). 
I/ES ”- 1 


( 2 . 2 ) 


It follows that 


1 


= inf 

xES n ~ 


f*(x), 


-= sup f*(x). 

xes "- 1 


One easily shows that /(u) = sup{a; • v : x £ K} and /*(t) = inf{A : ? € K}. This also implies that 
B ms C K C Bm # , and so if \K\ = 1, then m$\B\ < 1 < M£\B\. As mentioned in the introduction, 
the surface tension / and gauge function /* are dual in the sense that they satisfy a Cauchy-Schwarz-type 
inequality, called the Fenchel inequality: 

x ■ v < f*(x)f(v) 


for all x, v £ M". We may characterize the equality cases in the Fenchel inequality: for any u, x-v = /*(s)/(u) 
if and only if v is normal to a supporting hyperplane of K at the point £ dK. Indeed, v is normal to 
a supporting hyperplane of K at x £ dK if and only if v ■ (y — x) < 0 (so v ■ y < v ■ x) for all y £ K. This 
holds if and only if u ■ x = supjj/ • v : y £ K } = /(u). In particular, if re £ d*K, then /*( x) = 1 and 


f(v K (x)) = x ■ v K {x). 


(2.3) 
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We may compute the gradient of /* at points of differentiability using the Fenchel inequality. The gauge 
function /* is differentiable at Xg £ R™ if there is a unique supporting hyperplane to K at € dK. 

For such an x 0 > let vq = vk ( / ^) ) £ R™ be normal to the supporting hyperplane to K at ^ , so 

• vq = /(^o) by (2.3). We define the Fenchel deficit functional by G(x) = f{^o)f*{x) — x ■ vq. By the 
Fenchel inequality, G{x) > 0 for all x and G{x o) = 0, so G has a local minimum at xq and thus 

0 = VG(i 0 ) = /NV/,(*o)-m. 

Rearranging, we obtain V/*(xo) = jr^j- The 1-homogeneity of / then implies that 

/(V/*(®)) = 1. (2.4) 

Furthermore, this implies that 

x • V/*(x) = x ■ v k ( t = /*(z) (2.5) 

(alternatively, this follows from Euler’s identity for homogeneous functions). An analogous argument ensures 

that 

Vf(v K (x)) = x (2.6) 

for x £ d* K. Furthermore, we compute 

div7U = Wh (2. 7) 


Indeed, 


x trVx 

div —r = „ , s + x • V 


x n — 1 

f*(x) f*(x)' 

1 


f*(x) f*(x) 

where the final equality follows from (|2.5|). 


( f*(x)) 


V f*(x) n — 1 


f*(x) f*(x) 2 f*{x)‘ 


2.3. Properties of /3$, and 7 $. Using the divergence theorem, by approximation and the dominated 
convergence theorem, and (2.7), we find that for any y £ R n , 

dx 


x-y 


We may then write 
where 7 <s>(E) is defined by 


ld*E f*(x - y) 

P*(E)' 2 = 


v E {x)d'H n =(n—1) 


‘e f* (a; - y) ' 


2 $(£?)-(n-l) 7 *(£?) 


7 $(£ i ) = sup / 
t/effi" 


77,|Ta | 1 /”-| J E , | 1 /”' / 
dx 


( 2 . 8 ) 

(2.9) 


v&JeMx-vY 

The supremum in (2.9) is attained, though perhaps not uniquely. If y £ R n is a point such that 

ME) = f -A—, 

Je j*(x - y) 

then we call y a center of E, and we denote by y E a generic center of E. The Wulff shape K has unique 
center y E = 0. Indeed, take any y £ R”, y ^ 0, and recall that K = {/*(x) < 1}. Then 


dx 


dx 


dx 


dx 


Ik f*(x) Jk f*(x ~ y) Jk f*{x) J K + y f*{x) 

f dx 


r dx 

ll<\(K+y) f*(x) 


> 


I(K+y)\K f*(x) J 
A similar argument verifies that if \E\ = \K\. then 

7<s>(E) < 7 $ (K). 

Moreover, (n — 1)7 <s>(K) = §(K) = n\I<\. 


ldx — 


1 dx = 0 . 


K\(K+y) 


l(K+y)\K 


( 2 . 10 ) 






























A STRONG FORM OF THE QUANTITATIVE WULFF INEQUALITY 


9 


The following continuity properties of $ and 7 $ will be useful. 

Proposition 2.1. Suppose that {E 3 } is a sequence of sets converging in L l to a set E, and suppose that 
{f J } is a sequence of surface tensions converging locally uniformly to f, with corresponding surface energies 
{} and 

(1) The following lower semicontinuity property holds: 

$(£) < liminf <S>j{Ej). 

j->oo 


(2) The function 7 $ defined in (2.9) is Holder continuous with respect to L 1 convergence of sets with 
Holder exponent equal to \jv!. In particular, 

ME) - ^(F)\ < ^lEAF^ 
n — 1 

for any two sets of finite perimeter E, F C M n . Moreover, 

lim 7 *AEj) = 7 *(E). 

3 ~* 00 

Proof. Proof of (1): From the divergence theorem and the characterization f{v) = sup{x • v : f*(x) < 1}, 
one finds that the surface energy of a set E is the anisotropic total variation of its characteristic function 
Xe- 

*i(Ej) = TV fj ( XEj ) := sup { J div Tdx \ T G C'g(lR n ,R n ), /» (T) < l}. (2.11) 

Let T G C'g(K",K n ) be a vector field such that /*(T) < 1 for all x G M n . Then we have 

divTdx= lim / di vTdx = lim ||/* (T)||z,~(r™) / divSj dx < liminf $j(Ej), 

: 3^°°Je, Je , •J-*' 00 


IE 


where we take Sj = T/||/*(T)|| £ «,fign). Taking the supremum over {T G Cf{ 
the result. 


l ) : /*(T) < 1}, we obtain 


Proof of (2); By 

ME) ~M F ) < [ 

Je 


dx 


dx 


< 


dx 


Letting r be such that \rK\ = \EAF\ and recalling (2.10), we have 


dx 


Ieaf f*( x ~ Ve) 


< 


E f*(x-y E ) JfMx-Ve) Jeaf f*{x-y E y 

n\K\ 


f dx <I>(rK) 

/ TTF = M rK ) = - 7 - 

IrK J*\ x ) n- 1 


i\K\r 


-\EAF\ 


1/n 


( 2 . 12 ) 


n — 1 n — 1 

Thus 7 <s>(E) — 7 <s>(F) < '^Si\EAF\ 1 l n ' . The analogous argument holds for 7 $(F) — 7 $(E), implying the 
Holder continuity of 7 $. 

For the second equation, we note that if f J —> f locally uniformly, then /* —>• /* locally uniformly and 
M$ j —► M$. The triangle inequality gives 

|7 Zj( E j) ~ M e )I < M(Ej) - M E j)\ + \M E j) - l<s>(E)\. 

The second term goes to zero by the Holder continuity that we have just shown. To bound the first term, 
let y Ej be a center of Ej with respect to the surface energy <& 3 . If 7 $,. (E() > 7 $(Ej), then 


0 < M(Ej) -M E j) < f 

Je 


1 


1 


dx = 


1 


1 


Ej fl(x- y Ej ) f*{x~ VEj ) 

= J Rn X(E J+ y Sj )\B'(0) 


’Ej+y E fl{x) f*{x) 


dx 


IB t ( 0 ) fi{x) f*{x) 


dx. 
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For e > 0 fixed, the first integral goes to zero as j -A oo. For the second integral, we have 


1 


1 


IB e ( 0 ) fi(x) f*(x) 


dx < 


M$ i + il-f$ 


dx < Ce n 


lB e (0) 


Taking e —► 0, we conclude that 7 <s> {Ej) — 7 §(Ej) —>• 0 as j —>• 00 . The case where 7 <s>AEj) < 7 <s>(Ej) is 

analogous. □ 


Remark 2.2. With sequences as in the hypothesis of Proposition 2.1 above, /3$ has the following lower 
semicontinuity property: 

P<s>{E) < liminf/3$ (Eh). 

j->oo 


This follows immediately from parts (1) and (2) of Proposition 2.1 and the decomposition in (2.8). 

Lemma 2.3. For every e > 0, there exists 77 > 0 such that if \FAK\ < r], then \y F \ < e for any center yF 
ofF. 


Proof. Suppose \KAFj —> 0. By the triangle inequality, 


/ K 


dx 

f*(x) 


< 


dx 


dx 


Ik f*(x) Jf, f*(x - y F] ) 


dx 


dx 


By (2.12), the first two terms on the right hand side go to zero as j —> 00 , implying that 

f dx 


Fi f* - VF 3 ) J K f* {x - yF, ) 
0 a 
dx 


dx 

Ik /*(x -Vf,)' 


Ji< f* O) Jk f* (x - Vf, )' 

Because K has unique center y F = 0, we conclude that y Fj \ —> 0. 


□ 


We now introduce the relative surface energy and the anisotropic coarea formula. Given an open set A 
and a set of finite perimeter E, the anisotropic surface energy of E relative to A is defined by 


<S>(E-,A)= J f{v E (x))dU n -\x). 


8*EnA 

For a Lipschitz function u : R n —> R and an open set E, the anisotropic coarea formula states that 


/ f(—’Vu(x))dx= f &({u > r}; E) dr. 
J E J 0 


The anisotropic coarea formula is proved in the same way as the coarea formula (see, for instance, m 
Theorem 13.1]), replacing the Euclidean norm with / and /* and using (2.11). When u is bounded by a 
constant C on E, then applying the anisotropic coarea formula to w = C — u yields 

[ f(S7u(x))dx= f f(—S7w(x))dx= I $({G — u> r}\E)dt 
J E J E JO 

= / $({■« < C — r}; E) dr = / $({zi <r};E)dr 

Jo Jo 

Moreover, approximating by simple functions, we may produce a weighted version: 

f f(Vu(x))g(f*(x))dx= I $({u < r};E)g(r)dr 
J E JO 

whenever g : R —> [0, 00 ] is a Borel function. We w ill frequently use this weighted version with u(x) = f*(x), 
E a bounded set, and g(r) = which, using (2.4), gives 


IE 


dx 

f*(x) 


*({f*{x)<ry,E) dr _ r°° 


< 1 >(rK ; E) 


dr. 


(2.13) 
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We conclude this section with the following Poincare-type inequality, which shows that (3<s>{E) controls 
a^(E) for all sets E of finite perimeter. 

Proposition 2.4. There exists a constant C(n) such that if E is a set of finite perimeter with 0 < \E\ < oo, 
then 

a$(A) + S^E) 1 ' 2 < C(n)P*(E). (2.14) 

Proof. We follow the proof of the analogous result for the perimeter in [22] ■ Due to the scaling and translation 
invariance of a<j>,/3$, and <5$, we may assume that \E\ = \K\ = 1 and that E has center zero. We have 


mk)-me) = f tt-) - f rh= l ttt- / 

Jk f*\ x ) Je f *\ x ) JK\E f*\ x ) J E 
Therefore, adding and subtracting 4>(A")/n = (n — l)j$(K)/n in (|2.8|), we have 


dx 

Ie\k f*i x ) 


/3$(A) 2 = 6$(E) 


n — 1 


I' dx 

Ik\e f *( x ) 


f dx 
Ie\k f*i x ) 


We want to bound the final two integrals from below by a^(E) 2 . To this end, we let a := \E \ K\ = \K \ E\ 
and define the Tv-annuli Ari = Kr \ K and A l r = K \ K r , where R > 1 > r are chosen such that 
IArpI = \A\^ r \ = a. In particular, R = (1 + a) 1 /" and r = (1 — a) 1,/n . By (2.10) and (2.13), 

f 1 $(sT£T) 


dx 


> 


dx 


K\E f*( x ) JA 1iT f*i x ) 


and 


t E\K 


dx 

f*i x ) 


< 


dx 

' ! f *( X ) 


$(sA') 


ds = 


ds = 


ji i —2 


ds = 


n — 1 


[i- 


n— li 


ns n 2 ds = 


— [iZ"" 1 -!]. 


Subtracting the second from the first, we have 


n— 1 


f dx 
Ik\e f *( x ) 


f dx 
Ie\k f*( x ) 


> 2 - r n ~ l - R n ~ 1 . 


The function g(t) = (1 + t) 1 ^ n is function is strictly concave, with \(g(t) + g(s )) < g{\ + f) — C\t — s| 2 , 
and therefore 2 — [(1 + a) 1 /” + (1 — a) 1 /" ] > 8C|a| 2 . Thus 

/3<f(A) 2 > S^(E) + [2 - (1 - a) 1/n ’ - (1 + a) 1/n '} > 6*(E) + 8C|a| 2 = 5$ + 2C (\E \ K\ + | K \ E \) 2 

= 5g>(E) + 2C|A'AA| 2 > 6$(E) + 2Ca< s> {E) 2 . 

□ 


3. Stability for General Anisotropic Surface Energy $ 

In this section, we prove Theorem |1.1| We begin by introducing a few lemmas that are needed the proof. 
The first allows us to reduce the problem to sets contained in some fixed ball. 

Lemma 3.1. There exist constants Rq > 0 and C > 0 depending only on n and A T$ such that, given a set 
of finite perimeter E with \E\ = \K\, we may find a set E' such that \E'\ = \K\, E' C Br 0 , and 

ME) 2 < ME') 2 + C6<s>(E), 5$(E r ) < CS^,{E). (3.1) 

Proof. A simple adaptation of the proof of [29], Theorem 4.1] ensures that we may find constants So, Co, C\, 
and R 0 depending on n and such that C 0 S 0 <1/2 and the following holds: if <5$ (A) < <5 0 , then there 
exists a set E C E such that E C Bp o and 

\E\ > \K\(1 - Ci6$(E)), $(E) < $(A) + CoS^{E)\E\ l / n ' . 


(3.2) 
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If 6&(E) > Sq, then 


P%{E) < = 8${E) + 1 < ^-^S^E). 


n\K\ <5 0 

Simply taking E' = K, we have 6$(E') < 6$(E) and (3<s>(E) 2 < 7^(5$ (A), proving ( |3.1| ). 


On the other hand, if 5$(E) < So, let E' = rE with r > 1 such that \E'\ = \rE\ = \E\. By (2.8 1 , 

?/ , 2 ^(E)-^(E') n- 1 


P *{E) 2 -fe{E'y = 


n\K\ + n\K\ 
n — 1 


(ME')-ME)) 


< S*(E) + (r"-^*^) - ME)) ■ 


(3.3) 


Since E C E, 7 $(E I ) < 7 $(E), which implies that 

n— 1 


»i^i 




< ~TVT7-(r n_1 — l)7$(i?). 


By (2.10) and the fact that 7 $(K) = n\K\/(n — 1), 
n — 1 


l|/v| 


( r " _1 - 1)7* (£) < 


n — 1 
n|/v | 


»l^l 


(r"- 1 -l) 7 $(^)=r"- 1 -l, 


and since r > 1 , 


„n—1 


- 1 < r" - 1 = 


The first part of (3.2) implies that 


\E\-\E\ < Cx&*(E) < 


1^1 

Cr 


We have therefore shown that 


n — 1 
n|A'| 


|£| - l-Ci^(A) « l-Cido 

Cr 


(r"" 1 7*(B)-7*(^)) 


< 


1 - Ci5 { 


6*(E). 

<5$ (A); 


this together with (3.3) concludes the proof of the first claim in (3.1). 


In the direction of the second claim in (3.1), the first and second parts of (3.2) respectively imply that 

HE) $(£’) + CoS^{E)\E\ 1 / n ' 


®(E') = r n ~ l H E ) < 


< 


{i-cME)y/"' ~ (i ~cME)y /n ' ' 

A Taylor expansion in <5$ (E) of the right hand side shows that 

$(£') < $(£?) + CoS^(E)\E\^ n ' + —Icifa (£)$(£) + 0{S^(E) 2 ) 

n 

< <f>(E) + C5®(E)<f>{E) 

for Jo chosen sufficiently small. Thus 

r *(&) - *(K) ^ HE') - HE) ^ CHE)6*(E) 

n\K\ n\K\ n\K\ 

since &(E) < <1 >(K) + n\K\So- Finally, since E C B^ and E' = rE with r < 1/(1 — C\5o 7”, we have 
E' C B Ro for Rq = rR 0 . 

Let us now consider the functional 


□ 


Q(E) = HE) 


\K\mq, 

8M$ 


\ME) 2 -e 2 \+A\\E\-\K\\ 


(3.4) 


with 0 < e < 1 and A > 0. 
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Lemma 3.2. A minimizer exists for the problem 

min {Q(E) : E C B Ro } 

for A > 4 n and e > 0 sufficiently small. Moreover, any minimizer F satisfies 

\K\ 


1*1 > 


2 ’ 


<f>(F) < 2n\K\. 


(3.5) 


Proof. Let Q = inf{Q(F) : E C B Ro }, and let {Fj} be a sequence such that Q(Fj) —► Q. Since Fj C B Ro 
and $(Fj) < 2 Q for j large enough, up to a subsequence, Fj —y F in L 1 for some F C B Ro . The lower 
semicontinuity of $ (Proposition |2.1[1)) ensures that $(F) < oo. 


We first show that |F| > 7 ^. For any ?y > 0, Q(Fj) < Q + \K\rj for j sufficiently large. Furthermore, 


Q < Q(K) = $(K) 


e 2 \K\m<& 
8 M$ 5 


SO 


w - |K| i * 1 (*<*> - ™ + c ^t) ^ 

for e and 77 sufficiently small. Therefore \Fj\ > implying that |F| > ^ as well. 

We now show that liminf Q{Fj ) > Q(F) , so F is a minimizer. Recalling (2.8), we have 

I TC I m .t. dl 

Q{Fj) = <S>{Fj) + 

> HFj) + 


|A'|m$ 

$(F j )-(77-l) 7 $(F J ) 2 

+ A \Fj\ 

-\k\\ 

8M<$> 

n\K\ 1 / n \Fj\ 1 / 71 ' 

|A'|m$ 

®{F) — (n — 1 ) 7 $ (F,) 2 

| AT] 771 $ 

<K* 1 ) - 

8M<$> 

77|A'| 1 /”|Fj| 1 /"' 

8M<$> 

n|A'| 1 /"|F ; -| 1 /"' 


+ A| |*11 — |A'| |. 

Let a = liminf i>(Fj). Up to a subsequence, we may take this limit inhmum to be a limit. By the lower 
semicontinuity of <E>, a > 4>(F). Furthermore, 7 $ is continuous by Proposition |2.1| (2), so 

“ Q(Fj) > Q(F) + (a - *(F)) - l« - *(F)| 

= Q(F) + («- $ (F))(l-JU^ 

2 1 /"' 


1 - 


771$ 


8t7Mj> 


> Q(*). 


> Q(F) + (a - 4>(F)) (l 
Finally, e < 1 and therefore $(F) < Q(F) < Q(K) < 2n\K\. 

The following lemma shows that a minimizer of 13.41) satisfies uniform density estimates. 


□ 


Lemma 3.3. Suppose F is a minimizer of Q(E) as defined in (3.f) among all sets E C B Ro . Then there 
exist ro > 0 depending on n , A, and \K\ and 0 < Co < 1/2 depending on n and A such that for any x £ d*F 
and for any r < r 0 , 


CQTTlg 

MS 


u) n r n < | B r (x) n F| < ( 1 — 


comg 

MS 


uj n r 


(3.6) 


Proof. We follow the standard argument for proving uniform density estimates for minimizers of perimeter 
functionals; see, for example, m Theorem 16.14]. The only difficulty arises when handling the term 
^m^-\/3<s>{E) 2 — e 2 \ in Q(E), as it scales like the surface energy. 

For any x 0 £ d*F, let r < r 0 , where r 0 is to be chosen later in the proof and r is chosen such that 


Ti n ~ 1 (d*F n dB r (x 0 )) = 0. 


(3.7) 
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This holds for almost every r > 0. Note that if (3.6 1 holds for almost every r < ro, then it must hold for 
all r < r 0 by continuity; it is therefore enough to consider r such that (3.7) holds. Let G = F\B r (x 0 ). For 
simplicity, we will use the notation B r for B r (xo)- Because F minimizes Q , 


$(F) 


\K\m$ 


/3$(F) 2 — e 2 | + A||F| — |/v|| < $(G) + 


\K\m$ 

8M$ 


/3$(G ) 2 — e 2 | + A||G| — |F|| 


and so rearranging and using the triangle inequality, we have 

\K\m<& 


<F(F) < <F(G) + 


8M<j 


|/3$(F) 2 — /3$(G) 2 | + A|Fn B r \. 


We subtract $(F;K” \ B r ) from both sides; this is the portion of the surface energy where d*F and d*G 
agree. We obtain 

Q(F-,B r )< f f(yB T .)d'H n ~ 1 + |/? 4 >(F ) 2 — /?$(G) 2 | + A|F n B r \. (3.8) 

JdB r nF 4 > 


Indeed, this holds because (3.7) implies that 


$(G) = $(F;R”\F r )+ f f(u Br ) 

JdBrnF 


dU n 


We must control the term ^1™"’ \(3$(F) 2 — /3$(G) 2 | and require a sharper bound than the one obtained 
using Holder continuity of 7 $ shown in Proposition |2.1[ 2). Indeed, we must show that the only contributions 
of this term are perimeter terms that match those in ( |3.8| ) and terms that scale like the volume and thus 
behave as higher order perturbations. We have 


|/3$(F) 2 — /3$(G) 2 | = 


1 


< 


n\K\ l / n 
2 $(F) 
n\K I 1 / 71 


$(F) - (n - 1 ) 7 $(F) $(G) - (n - 1 ) 7 $(G) 


\ F \l/n' 

\F\~ 1/n ' - IGI-V^I 


|G|V"' 

|$(F) - H>(G)| + (n - 1) | 7 $ (f) ~ 7 $(G)| 
n\K\ 1 / n \G\ 1 / n ' 


The function v(z) = 1 — (1 — z ) 1//n is convex and increasing with u(l) = 1 , hence v(z) < z for z € [ 0 , 1 ]. 
Thus, as |G| = |F| — |F fl B r \, 


11*1 


-l/ n ' _ | G |~l/n'j _ IGfj-l/n' _ \ F 


1 /n^ 


< 


|FnF r | 
|G| 1 /"'|F| ‘ 


(3.9) 


Since 2|F| > \K\ by (3.5), 4|G| > \K\ for r 0 sufficiently small depending on n, so the right hand side of (3.9) 
is bounded by 8\K\- l ~ 1 / n |F D B r \. The coefficient is bounded by 4|/\ j 1 /" thanks to (3.5), so 


2$(F) 

n \K\ 1 / n 


I* 1 ! 


— l/n |^f| — l/n 


< 32|F| _1 |F D B r \. 


(3.10) 


Therefore, by (3.10) and again using the facts that 4|G| > \K[, 2|F| > \K\, and m$/M$ < 1, we have shown 
that 


| K\m$ 

8Af<j, 


\P<s>(F) 2 — /3$(G) 2 | < 4|F n B r \ + 


|$(F) — $(G)| n — 1 

2 n M<f, 2 n 


| 7 *(F)- 7 *(G)|. 


For the term |$(F) — <b(G)|, using (3.7), we have 


|*(F)-$(G)| = 


f f{v F )d'H n ~ 1 — f f{v G )dU n - 1 <$(F;F r ) + [ f{ VBr )dU n ~\ 

J d*F J d*G JdB r nF 


(3.11) 


(3.12) 
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using (3.7) and the fact that d*F and d*G agree outside of B r . Similarly, for the term | 7 $(F) — 7 $(G)|, 
when 7 $(A) > 7 $(G), thanks to (3.7) we have 


7$(F) - 7$(G) < 


(x~y F ) -Vf(x) jnjn-l 


3*f f*(x-y F ) 


dH n ~ l - 


( X - y F ) ■ vg(x) 
3*G f*(x-y F ) 


dH n 


< — ($(F;B r ) 




>dB r nF 


f{vB r )dU n 


The analogous inequality holds when 7 $(G) > 7 $(F), so 

|7$(^)-7$(G)| < — U{F-,B r ) + [ f(vB r )dU n - 
m <S> \ JdB r n F 

Combining (3.11), ( 3.12[ ), and ( |3. 13 1, we have shown 

|/?$(A) 2 -/3$(G) 2 | <A\FnB r \ + U^(F-B r )+ [ f{v Br ) dH r 

Z \ JdB r DF 


\K\m<s> 

8117$ 


(3.13) 


(3.14) 


Combining (3.8) and (3.14) and rearranging, we have 


B r 


i$(F; B r ) <\f f{v Br ) dU n ~ x + (4 + A) \F n . 

z 1 JdB r nF 

Proceeding in the standard way, we add the term \ f gg nF f{v Br ) dTL n ~ 1 to both sides, which gives 
1${F n B r ) < 2 [ f(u Br ) dn n ~ x + (4 + A) \F n Br I . 

Z JdB r f\F 

By the Wulff inequality, $(F D B r ) > n\K\ l / n \F D B r I 1 /”', and for ro small enough depending on n, A, and 
\K\, we may absorb the last term on the right hand side to obtain 

n\K\ x / n \F n B r \ 1 / n ‘ 


< 2 


/ (y Br ) dT~i 


n— 1 


(3.15) 


!dB r nF 


Let u{r) = \FC\B r \, and thus u'(r) = TL n 1 (dB r CiF), so the right hand side above is bounded by 2 M$u'(r). 
Furthermore, lATj 1 /™ > m$, so (3.15) yields the differential inequality 


nm$ 

8M$ 


< vf(r)u(r) 


— l/n' _ 


= n(u 1 /n y. 


Integrating these quantities over the interval [0, r], we get 


rricpr 


< u{rf/ n = \B r nF\^ r ' 


8 M $ 

and taking the power n of both sides yields the lower density estimate. The upper density estimate is 
obtained by applying an analogous argument, using G = F U B r (xo) as a comparison set for Xq £ d*F and 
r < ro satisfying (3.7). □ 


The following lemma is a classical argument showing that a set that is close to K in L 1 and satisfies 
uniform density estimates is close to K in an L°° sense. 


Lemma 3.4. Suppose that F satisfies uniform density estimates as in (3.6). Then there exists C depending 
on m$/M$, n, and A such that 

hd (dF,dI\) n < C\FAK\, 


where hcl(-,-) is the Hausdorfif distance between sets. In particular, for any y > 0, there exists e > 0 such 
that if \FAK\ < e, then A'l-r; C F C Ki +V , where K a = aK. 
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Proof. Let d = hd (dF, dK). Then there is some x £ dF such that either B d (x) is contained entirely in the 
complement of K or B d (x) is entirely contained in K. If the first holds, then the lower density estimate in 
(3.6 1 implies that 

\FAK\>\FnB d (x)\> C -^$d n , 


while if the second holds, then the upper density estimate in (3.61 implies that 

\FAK\>\B d (x)\F\>^d n . 


□ 


We will make use of the following form of the Wulff inequality without a volume constraint. 

Lemma 3.5. Let Rq > diam(K) and A > n. Up to translation, the Wulff shape K is the unique minimizer 
of the functional 

$(E) + A||E| - |if11 

among all sets F C B r 0 . 

Proof. Let £ be a minimizer of <E >(F) + A||E| — |A'|| among all sets of finite perimeter F C -Br 0 ; this 
functional is lower semicontinuous so such a set exists. Comparing with K , we find that 

$(A) + A||A| - |if11 < $(if) =n|if|. (3.16) 

The Wulff inequality implies that \E\ < |if|, and so &(E) > ?r|E| 1 / n l/vl 1 /" > n\E\. Thus (3.16) implies 
that A (|if | — |2?|) < n(|if| — |E|). Since A > n, it follows that \E\ = \K\. It follows that E must be a 
translation of K , the unique (up to translation) equality case in the Wulff inequality. □ 

We are now ready to prove Theorem |1.1| 


Proof of Theorem 0 By ( |1.4| ), we need only to show that there exists a constant C = C(n) such that 

/3<t,(A;) 4n/(n+1) < C6*(E), (3.17) 

for any set of finite perimeter E with 0 < |2£| < oo. By Lemma [3. 1| it suffices to consider sets contained in 
Br 0 . Let us introduce the set 

-Mcj> 


?K = {f : 


for Af > 1, recalling M d , and defined in (2.2). In Steps 1-4, we prove that, for every Af > 1, there exists 
a constant C = C(n,Af) such that (3.17) holds for any surface energy $ corresponding to a surface tension 
/ G Wy. In Step 5, we remove the dependence of the constant on Af. 


Step 1: Set-up. 

Suppose for the sake of contradiction that (3.17) is false for some Af. We may then find a sequence of sets 
{Ej} with Ej C Br 0 and a sequence of surface energies {'I» y }, each with corresponding surface tension 
€ Fjq-, Wulff shape Kj, and support function /*, such that the following holds: 


l^-l = \Kj\ = 1, 

->• 0 , 

QjiEj) < &j(Kj) + ci/3c S>j (Ej) in ^ n+X \ (3.18) 

where c\ = C\ (Af, n ) is a constant to be chosen later in the proof. 

Each f J is in fFjg and is normalized to make Kj = 1 implying that {f 3 } is locally uniformly bounded 
above, and hence, by convexity, locally uniformly Lipschitz. By the Arzela-Ascoli theorem, up to a sub¬ 
sequence, f J -A /°° locally uniformly. The uniform convergence ensures that this limit function /°° is a 
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surface tension in T\j. We denote the corresponding surface energy by $ 00 , Wulff shape by K^ and support 
function by fffi. Note that (A^l = 1 . 

There exists c(Af) such that &j(E) > c(J\f)P(E) for any set of finite perimeter E , again thanks to 
fi G Fjg and \Kj\ = 1. Then, since $j(Ej) —> n (as $>j(Kj) = n ), the perimeters are uniformly bounded. 
Furthermore, EjC B Rq , so up to a subsequence, Ej —> E^ in L 1 with lA^J = 1. 


Proposition 2.1 T) implies that < 3? 00 (A , 00 ) < lim §j(Ej) = n, so by the Wulff inequality, E^ = K up 

’2) then ensures that lim r y$ j (Ej) = 7 $ oo (A' 00 ) = ^ 37 , and 


to translation. Furthermore, Proposition 
therefore, by ( 2 . 8 ), 


2.1 


lim Aj, (Ej) 2 = lim - ($j(Ej) - (n - 1 ) 7 $ . (Ej)) = 0. 

00 J i-ioo n 


Step 2: Replace each Ej with a minimizer Fj. 

As in pi], the idea is to replace each Ej with a set Fj for which we can say more about the regularity. We 
let €j = p^^Ej) and let Fj be a minimizer to the problem 


min Qj(F) = <f>j(F) + 


8M* 


\fcSF) 2 -e 2 A+A\\F\-l\ : FcB 


Rq 


for a hxed A > 4n. Lemma 3.2 ensures that such a minimizer exists. As before, <&j(Fj) > c{Af)P{Fj). 
Pairing this with (3.5) provides a uniform bound on P(Fj). so by compactness, Fj —> F 0 0 in L 1 up to a 
subsequence for some F x C B Rq . 

For each j, we use the fact that Fj minimizes Qj , choosing Ej as a comparison set. This, combined with 
(3.18) and Lemma 3.5 yields 


1 


+ W \P*A F :Y - S i + A |l*il - 1 | < Qi(Fj) < *j&i) 

< $j{Kj) +Cl ef /(n+1) < $j{Fj) + A||Fj| - l| + Cl e 4n/(n+1) . 


(3.19) 


It follows that ^f\P'S> j {Fj) 2 — e|| < cie^ n ^ n+1 \ immediately implying that P^, j (Fj) —> 0. Moreover, rear¬ 
ranging and using the fact that e ? 0 and > 2 , we have 

SlJvST < </MU> 2 , 

where the exponent (n + l)/2n is chosen so that, taking the power 2 n/(n + 1 ), we obtain 

4n/(n+i) < 2 fe.{Fj) 4n ^ n+1 \ (3.20) 


In the last inequality in (3.19), if we replace Fj with arbitrary set of finite perimeter E C B Ro , then we 
obtain 

<f>j(Fj) + A\\Fj\ -1| < $,-(£?)+A||£7| - l| + Cl e 4n/(n+1) , 

again using Lemma |3.5| Taking the limit inferior as j —► 00 , this implies that F a0 is a minimizer of the 
problem 

min{$ 00 (F)+A||A|-l| : F C B Ro } , 


and so F^ = K a Q up to a translation by Lemma 3.5 With no loss of generality, we translate each Fj such 
that inf{ | (Fj + z)AK oa |:z 6 M"} = | FjAK x \. 

Step 3: For j sufficiently large, \Kj C Fj C 2 I\j and \Fj\ = 1. 

Lemma H implies that each Fj satisfies uniform density estimates, and thus for j sufficiently large, 
Lemma 3.4 ensures that \Kj C Fj C 2 Kj, as \KjAFj\ < \KjAKoo] + AAc AA ? and both terms on the right 


hand side go to zero. 
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3.2 


Let r : j > 0 be such that \rjFj\ = 1 . We may take r 3 Fj as a comparison set for Fj ; r ; - < 2 by Lemma 
so VjFj C 4Kj C B Ro as long as R 0 > 4> CJ\f, the second inequality following from \Kj\ = 1. Since 
/3<i> is invariant under scaling, Qj(Fj ) < Qj(rjFj) yields 

* j {F j )+A\l-\F j \\<r'r 1 * j (F j ). 


This immediately implies that r :j > 1 for all j, in other words, \Fj\ < 1. Furthermore, Tj 


(3.21) 
1 because 


Fj 


Kra in L 1 and I K a 


3 \ — / j 

= 1. Suppose that, for some subsequence, r,- > 1. Then, using \Fj = 1 / r" , (|3.21) 


implies 


A < 


W 1 - 1) 


r J~ 1 


<W')- 


(3.22) 


For any 0 < 77 < ^ and for j sufficiently large, the right hand side is bounded by (1 — rj)$j(Fj), as 
lim 

r —>■!+ 


r"(r" -1) _ 7l-l 


_ n i^ F ur thermore, <frj(Fj) < n + e 2 since Qj(Fj) < Qj(Kj), so (3.22) implies that 


A < (1 — r])$j(Fj) < (1 - rj) (n + e 2 ) < n 

for j sufficiently large. Since n < A, we reach a contradiction, concluding that \Fj\ = 1 for j sufficiently large. 
Step 4 : Derive a contradiction to ( |3T8| . 

We will show that /3^, j (Fj)' ln ^ n+1 ^ < (Fj ), which in turn will be used to contradict (3.18). Adding and 

subtracting the term $j(Kj)/n = (n — 1 ) 74 ^ (Kj)/n to (2.8), we have 


A F :) 2 < 




— ( Fj ) + 


n — 1 
n 

n — 1 
n 


f dx 

If i f*(x) 


— (Fj) + 


n — 1 


da; 


'fAK, 


1 -t- dx 

f*(x) 


Ikj fl(x) 

1 


r dx 
If j fZ(x) 


Kj\Fj fl(x) 


— 1 dx 


We now control the last term in terms of 6^> j (Fj). Note the following: since \Kj C Fj C 2ATj, the last term 
above is bounded by C\FjAKj\ < S^^Fj) 1 ! 2 . This could establish (3.17) with the exponent 4. However, 
with the following argument, we obtain the improved exponent 4 n/(n + 1). 

As noted before, Lemma ( 


3.3 


1 implies that each Fj satisfies uniform density estimates (3.6 1 with /M<^ j > 
l/N. The lower density estimate provides information about how far /* (a:) can deviate from 1 for x £ Fj\Kj, 
thus bounding the first integrand. Indeed, arguing as in the proof of Lemma 3.4 for any x £ Fj \ Kj , let 
d= fi (x) — 1. The intersection Kj n B d (x) is empty by the definition of /* , and thus Fj n B d (x) C Fj \ Kj . 
Therefore, for x £ d*Fj \ Kj 


31 


^d n < \B d (x) D Fj\ < \FjAKj\ < (Fj) 1 ! 2 


by the lower density estimate in (3.6) and the quantitative Wulff inequality as in (1.4). In fact, this bound 
holds for any x £ Fj \ Kj ; since Fj is bounded, for any x £ Fj\ Kj , there is some y £ d*Fj \ Kj such that 
/*( x) < fUy). Therefore, fi(x) — 1 < C<5$. (Fj) 1 ! 271 for all x £ Fj\ Kj, and so 


1 - 


1 


<Fj\Kj 


f* (x) 


dx < f f*(x) — 1 dx < f C5<p.(Fj) ll,2n dx 
JFj\K, J fak , 


Fj\Kj 


(3.23) 


= C\FjAKj\5^(FjY/ 2n < C6^(Fjf/ 2 +^ 


where C = C(J\f,n) and the final inequality uses ( |1.4| ) once more. The analogous argument using the upper 
density estimate in (3.6), paired with the fact that eventually 1 Kj C Fj, provides an upper bound for the 


size of 1 — /* (x) for x £ Kj \ Fj, giving 

f 1 


' K i\Fj fl(x) 


— ldx<2 


[ 1 - fi (x) dx < C8 *. (ii» 1/2+1/2n . 

JKj\Fi 


(3.24) 
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Combining (3.23) and (3.24), we conclude that 


/M^-) 4n/(n+1) < Cm hAFj) 


(3.25) 


where C\ = 


We now use the minimality of Fj. comparing against Ej, along with (3.18) and (3.20) to obtain 

+c ie f /( " +1) < + 2ci/3<s>AFj) 4n /( n+1 \ 


By (3.20), /3 <i, j (Fj) is positive, so by choosing c\ < n/2C \, this contradicts (3.25), thus proving (3.17) for 
the class fFj\f with the constant C depending on n and A f. 


Step 5: Remove the dependence on Af of the constant in (3.11). 

We argue as in [15]. We will use the following notation: is the surface energy with Wulff shape K , surface 

tension f K , and support function fff. We use 5k, Pk , and 7 k to denote 5<s, K ,p^, K , and respectively. 

By John’s Lemma ([28) Theorem III]), for any convex set K C R n , there exists an affine transformation 
L such that detL > 0 and B\ C L{I\) C B n . This implies that M L ( K )/ m L(K) < n and so f L ^ K > £ F n . Our 
goal is therefore to show that /3$(13) and 5$(E) are invariant under affine transformations. Indeed, once we 
verify that Pk{E) = p L{K ){L(E)) and 5 K {E) = 6 L ( K ){L(E)), we have 

Pk(E )‘ 1 " /( ” +1) = p i( K)(L(i3)) 4n/(n+1) < C(n)6 L{K) (L(E)) = C(n)5 K (E), 
and (3.17) is proven with a constant depending only on n. 


Suppose E is a smooth, open, bounded set. Then 


$k(E) = lim 

e ->0 


\E + cK\ - \E\ 


this is shown by applying the anisotropic coarea formula to the function 

inf{/*(a: — y) : y £ dE} if x £ E c 
- inf{/*(x — y)'.y£ dE} if x £ E 


d K (x, dE) := 


and noting that {E + eK ) \ E = {x : 0 < d K (x, dE) < e}. 

Since L is affine, | L(E + eK)\ — \L(E)\ = det L (| E + cK\ — 1131), and so 

ME) = lim + = 11 wItW) 

e->o e det L det L 

Since \E\ = |L(13)|/detL, we have 

®k{E) 1 _ ®l(k){L(E)) 


Sk ^ n\K\ 1 / n \E\ 1 l n ' 1 n\L{K)\ 1 / n \L{E)\ 1 / n ' 1 Sl F)( l ( e )), 


and thus 5k{E) is invariant. Similarly, 

L~ 1 (z) — y 


f*{L 1 z- y ) = inf j A : 
and thus 


A 


dx 


£ K} = inf | 

dz 


A 


z - L(y) 

A 


L{K))=f^ K) ( z-L(y )), 


dz 


Je f*{x - y) 

Taking the supremum over y £ 


/ l(e) f* ( L Hz) - y) det L J L(E) fF K ) ^ _ £( y )) deti 
i n of both sides, we have 

7 l(k){L(E)) 


7 k{E) = 


det L 


From (2.8), 


fe(£) = (IdaWiWisMV' 2 


n\K\ 1 / n \E\ 1 / n ' 
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We have just shown that, for the denominator, 


1 


n\K\Hn\E\Hn' 


1/2 


and for the numerator, 


$ K (E)-(n-l) lK (E) 


1/2 


/ det L \ 1 ^ 2 

= ' v n|L(A')| 1 /«|L(i;)|i/™' ) ’ 

<fr L(g) (£(£)) ^ (n - 1)7l(a) (L(E)) \ 1/2 
det L / 


The term det L cancels, yielding 

O in ( <S>L(K){L{E))-{n^ 1)7L(A) \ 1/2 

PkK } ^ n|L(Ji)|l/n|L( S )|l/n' J 

showing that /3k(E) too is invariant. 


Pl(k)(L(E)), 


□ 


4. The case of A-elliptic Surface Tension 


In this section, we prove Theorem 1.3 This proof closely follows the proof of (1.5) in 


Using a 

selection principle argument and the regularity theory for (A, ro)-minimizers of <E>, we reduce to the case of 
sets that are small C ' 1 perturbations of the Wulff shape K. In }22| . this argument brings Fusco and Julin to 
the case of nearly spherical sets, at which point they call upon (1.13), where Fuglede proved precisely this 
case in [ 21 ]. 

We therefore prove in Proposition 1.9 an analogue of (1.13) in the case of the anisotropic surface energy $ 
when / is a A-elliptic surface tension. The following lemma shows that if FI is a small C 1 perturbation of the 
Wulff shape K with | E\ = \K\. then the Taylor expansion of the surface energy vanishes at first order and 
takes the form £3- We then use the quantitative Wulff inequality as in (|1.4[) and the barycenter constraint 


along with (4.1) to prove Proposition 1.9 


Lemma 4.1. Suppose that $ is a surface energy corresponding to a X-elliptic surface tension f, and E is 
a set such that \E\ = \K\ and 

dE = {x + u(x)uk(x) : x £ dK} 

where u : dK —> R and HuHcqaA') = e - There exists eo > 0 depending on A and n such that if e < eg, 

$(E) = $(K)+l [ (Vu) T V 2 f(y K )VudH n - 1 -l f H K u 2 dU n ~ x + e 0(||«|& 1(Wf) ), (4.1) 

1 JdK 1 JdK 

where Hk is the mean curvature of K and all derivatives are restricted to the tangential directions. 


Remark 4.2. The second fundamental form IIk of K satisfies 

V 2 /(pr:( 2 ;))IIk(x) = Wt x ok for all x £ dK. 

Therefore, Hk = tr(IlK) is equal to tr(V 2 /IlK) and thus ( |4.1| ) agrees with, for example, [9] Corollary 4.2], 

Proof of Lemma \4~l | For a point x £ dK, let {ti, ..., t„_i} be normalized eigenvectors of Vpr-, where each 
Ti corresponds to the eigenvalue A^. This set is an orthonormal basis for T X K, and thus {n,..., r n _i, vk\ 
is an orthonormal basis for R n . A basis for T x+UVk E is given by the set {g\,... ,g n _i}, where, adopting the 
notation iq = d Ti u, 

gi = d Ti [x + upk] = (1 + A iU)Ti + UiP K - 

We make the standard identification of an (n — l)-vector with a vector in R" in the following way. The 
norm of an (n—l)-vector v\A- ■ is given by |uiA.. .Au ra _i| = | det(ui,... ,v n -i)\. If |iqA.. .Au ra _i| 0, 

then the vectors Vi,..., v n -\ are linearly independent and we may consider the n— 1 dimensional hyperplane 
II spanned by v\, ■ ■ ■ , v n -\. Letting v be a normal vector to II, we make the identification 

id A ... A Vn— l = ±|m A ... A v n - 1 | v. 
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where the sign is chosen such that det(ui,..., v n —i, ±u) > 0. In particular, we make the identifications 

ST A ... A g n -1 


n A ... A t„_ i = v K , 


= ue, and tl A ... A vk A ... A t„_i = —r». 


|<?i A ... Aflf„_i| 

The sign is negative in the third identification because 

det(Ti,..., v K ,..., t„_i, —Tj) = — det(ri,..., —t*,..., t„_i, v k ) = det(r 1( ..., t u ... v K ) = 1. 

We let w := g\ A ... A g n -i, and so 

w = [(1 + A iu)n + u\v K ] A ... A [(1 + X n -iu)T n -i + u n -iv K \ 

n —1 n —1 

= n (1 + A iu)v K ~^2 u i JJ(1 + A ju)Ti 

i=l ijtj 


i =1 


1 + H k u + ^ XiXjU 2 vk - ^2 Ui 1 + ^2 \j u T i + e 0(\ u \ 2 + l^^l 2 )- 

i<j i =1 jjii 


(4.2) 


In order to show (4.1|, the volume constraint is used to show that the first order terms in the Taylor 
expansion of the surface tension vanish. We achieve this by expanding the volume in two different ways. 
The divergence theorem implies that 

n\E\ = [ x-VEd'H n ~ 1 = f (x + uvk) ■ f 7 —r^) M = / (x + u vk) ■ w d'H n ~ 1 . 

JdE JdK \\ w \) JdK 

Adding and subtracting vk = t 1 A ... A r n _i, and using (4.2) and the fact that vk • ly = 0, we have 

n\E\= f x ■ v K d'H n ~ 1 + f u +x ■ (w - v K ) + H K u 2 dTi 71 - 1 + eO(\\u\\ 2 H i, dK )). 

JdK JdK 

Since J gK x ■ vk dH n ~ l = n\K\, the volume constraint |2£| = \K\ implies that 

f x ■ (w — v K ) d'H n ~ 1 = - [ u +H K u 2 d'H n ^ 1 + eO(\\u\\ 2 H i, gK) ). (4.3) 

JdK JdK 

Now we expand the volume in a different way. Because / is a A-elliptic surface tension, the Wulff shape 
K is C 2 with mean curvature depending on A and n. Therefore, there exists to = to{X 7 n) > 0 such that the 
neighborhood 

D = {x + tv K (x) : x S dK, t € (— t 0 ,t 0 )} 

satisfies the following property: for each y £ D, there is a unique projection 7r : D —> dK such that 7 r(y) = x 
if and only if y = x + tvx{x) for some t £ {—to, to)- In this way, we extend the normal vector field vk to a 
vector field N K defined on D by letting N K : D —> M n be defined by N K {y) = v K (i r(y)). We also extend u 
to be defined on D by letting u(y) = u(n(y)) for all y £ D. Therefore, if e 0 < to, dE may be realized as the 
time t = 1 image of dK under the flow defined by 

= uN K (i/)t(x)), ipo{x) = x. 

Such a flow is given by i^t(x) = x + tuNx , and so Vi/> t (a:) = Id + tA where A = X/(uNk)- A small adaptation 
of the proof of [30j Lemma 17.4] gives 


Jipt = 1 +t tr(A) + —(tr(A) 2 - tr(A 2 )) + eO(|u| 2 + |Vu| 2 ). 


(4.4) 
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Integrating by parts, it is easily verified that 

[ tr (A) 2 - tr (A 2 ) dx = [ div (■ uN K div (uN K )) dx - [ V {uN K )^di{uN K )^dU^ 1 
Jk Jk JdK ii=x 

= f di v(uNk div (uNk)) dx — f uS7u ■ uk dK n_1 . 

Jk JdK 

The second equality is clear by choosing the basis Ti,... , r„_i, r„, where r n = v E - Furthermore, the 
divergence theorem implies that 

I div (uNk div (uNk)) dx = j u div {uN k) dK n ~ x = j uVu- vk + Hku 2 d'H n ~ 1 , 

Jk JdK JdK 


so that 


Therefore, the volume constraint \K\ = \E\ = \ipi(K)\ implies that 

1 


[ tr(A ) 2 — tr(A 2 ) dx = ( H k u 2 d'H r ‘ 

Jk JdK 

With this and ( |4.4| ) in hand, we have the following expansion of the volume: 

\MK)\ = f J^tdx= \K\+t f ud'H n - 1 + t - f H K u 2 dH n - 1 +t 3 eO{\\u\\ 2 H1 , dK) ) 
Jk JdK 2 JdK 

H K v?dH n - l + eO(\\u\\ 2 H i iBK) ). 


/ udK n ~ l =-- 

IdK 2 


'dK 


jn— 1 


Combining ( |4.3[ ) and ( |4.5| ), we conclude that 

[ x ■ (w - v K ) dU n ~ x = [ u 2 H K d'H n ~ 1 + eO(\\u\\ 2 H i, dK) ). 

JdK 2 J aK 

We now proceed with a Taylor expansion of the surface energy of E: 

*(E)=[ f(v E )dU n ~ l =( f(J^-)\w\d / H n ~ 1 = f f(w) dK r 

Jd-E JdK ''M' JdK 

= f f{y K ) dU^ + [ Vf(v K ) ■ (w - v K ) dU n ~ x 

JdK JdK 

+ \ [ [w-v K ] T V 2 f(v K )[iv-is K }d'H n - 1 + eO(\\u\\ 2 H1{dK) ), 

2 JdK 

so, recalling that V/(vk(x)) = x by ( |2.6| ), 

Th - 1 

®(E) = ${K) + [ x ■ (w - u K ) dW 1 " 1 + \ I u iUj{Tj^ 2 f{v K ) T j)d'H n - 1 +e°{\\u\\ 2 m{dK) ). 

JdK 2 j dK . , =i 


(4.5) 

(4.6) 


Applying (4.6 1 yields (4.1), completing the proof. 


□ 


We now prove Proposition 1.9 using (4.1) as a major tool. 


Proof of Proposition HB Suppose E is a set as in the hypothesis of the proposition, i.e., |£?| = \K |, bar A = 
bar A', and 

dE = {x + u(x)vk{%) ■ x G dK}, 

where u : dK —► R. is a function such that u G C 1 (dK) and ||u||c' i (eji’) = e < ei with ei to be fixed during 
the proof. Up to multiplying / by a constant, which changes A by the same factor and leaves 
unchanged, we may assume that \K\ = 1. Let 


B{u)=\ f (Vu) T S7 2 f{ VK )S7udU n ~' [ H k u 

2 JdK 2 J aK 


dK 


n— 1 
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so that, by (4.1) 


5®(E) — ~B{u) + e 0(\\u\\ 2 H i( aK }) 


(4.7) 


as long as e\ < eo for £o from Lemma [4.1| 

Step 1: There exists C = (7(n, A, to$/Mj>) such that, for e± small enough depending on m$/M$ and X, 


Uj 


\u\dH r 


< C6${E). 


(4.8) 


Step 1(a): There exists C = C{n, m$/M$) such that, for e\ = ei(m<s,/M$) small enough, 

\EAK\ < CS^E) 1 ' 2 . (4.9) 

The quantitative Wulff inequality in the form ( |1.4| ) states that \EA(K + Xo)| < C(n)5<s,(E) 1 / 2 for some 
Xq eR", so by the triangle inequality, 

\EAK\ < C’(n)<5$(£l) 1/2 + \{K + x 0 )AK\. (4.10) 

It therefore suffices to show that \(K + Xq)AK\ < CS^(E) 1 / 2 . By [30) Lemma 17.9], 


2n 

\I\A(K + x 0 )| < 2|x 0 | P(K) < -|x 0 |. 


(4.11) 


Furthermore, the barycenter constraint bar if = bar A' implies that 

x 0 = / Xq dx = / x dx — / x — Xq dx = / x dx — / x dx. 

Jk Je Jk Je Jk+x 0 

For ei small enough depending on E, K + Xq C a fact that is verified geometrically since 

|xq| —)■ 0 as e — y 0 and thus |rco | may be taken as small as needed. Therefore, 


l*o I = 


/ x dx — / xdx 

IE Jk+x o 


< 2M$\EA(K + x 0 )\ < M^Ci+d^E) 1 / 2 , 


where the second inequality comes from (1.4). This, (4.11), and (4.10) prove (4.9). 


Step 1(b): For e\ sufficiently small depending on A and n, 

[ \u\d'H n ~ 1 <2\EAK\. (4.12) 

JdK 

Let dx{x) = dist(x, dK). As in the proof of Lemma |4~I| there exists to = t 0 (X,n) such that for all t < to, 

{dx = t} = {x + tux(x)}. Take ei < to and let G t = {dx =t}n(E\ K). Then 

E\K = {x + tvx : dK : u(x) > 0}, t £ (0, m(x))}, 

G t = {x + tvx : x £ {x £ dK : u{x) > t}}. 

The coarea formula and the area formula imply that 

p poo p poo p 

\E\K\= \Vd K \dx= dt dU n ~ l = / dt J(Id + tv K ) dU n ~ l , 

Je\k Jo Jg ± Jo J{u>t} 


so 


mi 


dt 




' {u>£} 


If |{« > 0 \ dt = \ [ u + dU n K 

z JO 1 JdK 


The analogous argument yields |A'\A| > | f dK u dTL n 1 , and (4.12) is shown. Combining (4.9) and (4.12) 


implies (4.8). 
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Step 2: There exists C = C(n, \,m$/M$) such that, for e\ = ci(n, A, m<s>/M<s,) small enough, 

IMI ‘h i (ok) — C6$(u). 

The A-ellipticity of / implies 


(4.13) 


IdK 


\Vu\ 2 dH r ' 


~ A. 


[ (\7u) T V 2 f(u K )(Vu)dH n - 1 = \hB(u) + [ H K \u\ 2 <TH r 

JdK a V J dK 

The Wulff shape K is bounded and C 2 , so H E is bounded by a constant C = C(n, A). Therefore, 

[ \S7u\ 2 dU n - 1 < \b{u) + C f | u\ 2 dH n -\ (4.14) 

JdK a JdR- 

As pointed out in [12l proof of Theorem 4], one may produce a version of Nash’s inequality on dK that 
takes the form 

f \u\ 2 dH^ 1 <c^ n+2)/n [ \S7u\ 2 dU n - 1 + c J f \ u \dU n - 1 ) 2 (4.15) 

JdK JdK V {n+)/ K JdK J 

for all 77 > 0, where c is a constant depending on Hk (and therefore on A and n) and Indeed, as in 

R ra , this form of Nash’s inequality is a consequence of the Sobolev inequality on dK, shown in [33J Section 
18]. We pair (4.15) with (4.14) and (4.8) to obtain 

/ | Vu| 2 dU n ~ x < \b(u) + Cii n+2)/n f | Vu| 2 <M 
JdK A J dK 


n—1 


c 


fj( n + 2)/2 


S*(E). 


For r] small enough, we absorb the middle term into the left hand side. Then, recalling (4.7), we have 

1 


/ \\/u\ 2 dH n - x < CS^(E) + eO(\\u\\ 2 Hl(dK) ). 

JdK 


Combining this estimate with (4.15) and (4.8), we find that J gK \u\ 2 dTL™ 1 is also bounded by C6$(E) + 
e °{\\ u \\ 2 HHdK))- Therefore, 


Ml 2 h 1 (8K) — CS^(E) + e 0 (IMIffi(aic)) • 


Finally, taking e\ small enough, we absorb the second term on the right, proving (4.13). 


□ 


We now show that if dE = {x + uv K : x € dK} with IMIcqaif) small, then (3<&(E ) is controlled by 
M \H 1 (dK)- With the notation from the proof of Lemma 14. ll 


n\K\fe(E) 2 < [ f{v E ) - 

JdE 


f*(x) 


V E dH n - x = 


1 = / fW 

JdK 


— x ■ w dTL r 


From the expansion of 4> in the proof of Lemma 4.1 and the fact that x ■ v E = J{vk) by (2.3), the right 
hand side is equal to 

^4C (VU)TV2/( ^ )VW ^ n " 1+e0( l|u|l ^ 1( ^ )) - C ' l|u|l « 1( ^ )+e0(l|u|l « 1( ^ )) ’ 

where C = \\V 2 f\\c°(dK)- F° r e sufficiently small, we absorb the term e 0{\\u\\ 2 Hl ^ dK ^) and have 


P<s>{E) 2 < 


C 

n\K\ 


mm 


H'idK)' 


(4.16) 


Remark 4. 3. T his is the first point at which we use the upper bound on the Hessian of /. In other words, 
Proposition 1.9 still holds for surface tensions / £ C' 1 ’ 1 (K™\{0}) that satisfy the lower bound on the Hessian 
in the definition of A-ellipticity. 

Next, we prove Theorem |1.3[ for which we need the following definition. 
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Definition 4.4. A set of finite perimeter E is a (A,ro)-minimizer o/4>, for some 0 < A < oo and rg > 0, 

if 

<f>(E; B(x, r)) < <f>(F ; B(x, r)) + A|£AF| 

for EAF CC B(x, r) and r < r$. 


Proof of Theorem [77 


Proposition 2.4 implies that the proof reduces to showing 

/3 $(E) 2 < C6$(E). 


(4.17) 


where C = C(n, A, || V 2 f\\c°(dK), m$/M<p). Suppose for contradiction that (4.17) fails. There exists a 
sequence { Ej } such that \Ej\ = \K\ for all j, 5$(Ej) 0, and 


^(Ej) < $(A') + c 2 /3 $ (^) 2 


(4.18) 


for C 2 to be chosen at the end of this proof. Arguing as in the proof of Theorem |1.1| we determine that, up 
to a subsequence, {Ej} converges in L 1 to a translation of K. As in the proof of Theorem (and as in 
[22]), we replace the sequence {Ej} with a new sequence {Fj}, where each Fj is a minimizer of the problem 


min < Qj(E ) = &(E) + 


\K\m<& 

8 M$ 


ME)'- ej] + A\\E\ ~ \K\\ : E C B Ro 

Continuing as in the proof of 


with €j = P$(Ej)', existence for this problem is shown in Lemma 3.2 
Theorem 0 we determine that 

e 2 < 2/3 *(Fj) 2 , (4.19) 

that up to a subsequence and translation, Fj —► K in L 1 , and that \Fj\ = \K\ for j sufficiently large. By 


Lemma |3.3[ each Fj satisfies uniform density estimates, and so by Lemma |3.4[ for any 77 > 0, we may choose 
j sufficiently large such that Ki_ n C Fj C K\ +v . 

Arguing as in j22] . we show that Fj is a (A, ro)-minimizer of $ for j large enough, where A and r 0 are 
uniform in j. Let G such that GAFj CC B r {x 0 ) for Xq £ Fj and for r < ro, where ro is to be fixed during 
the proof. For any 77 > 0, if B r (x 0 ) C Ad-r,, then trivially 47(G) > d>(Fj). If B r (x 0 ) <£. then for 

77 sufficiently small, Lemma 2.3 implies that \yp. | <1/4 and \hg\ <1/4. Furthermore, by choosing 77 and 
ro sufficiently small, we may take B r {X q) n K 1 / 2 = 0. The minimality of Fj implies Q(Fj ) < Q(G ); after 
rearranging and applying the triangle inequality, this implies that 

|/L|t7T$ 


<S>(Fj)<<f>(G) + A\FjAG\ 


8M$ 


/3$(G) 2 - fo{Fj 


(4.20) 


As in (3.11) in the proof of Lemma 3.3 

\K\m<s> 


8 Ms, 


W $(F) 2 -^(G) 2 < 


l< $(Fj)-*{G)\ , 1 7 *(-^j') — 7$(G)| 


MFjAG\ 


for ro small enough depending on n. If 47 (A)) < 47(G), then the (A, ro)-minimizer condition is automatically 
satisfied. Otherwise, subtracting 14>(F)) from both sides of ( |4.20| ) and renormalizing, we have 

4>(A,) < 47 (G) + | 7 *(G) - MFj)\ + (8 + 2A)|A,AG|. (4.21) 

To control jy$(G) — 7 $(F))|, we need something sharper than the Holder modulus of continuity of 7 $ given 
in Proposition |2.1[ 2). Indeed, 7 $ is Lipschitz continuous for sets whose intersection contains a ball around 
their centers: 


74>(-P/)-74>(G) < [ 
J F 


dx 


dx 


dx 


F, /* 0 - VFj ) J G f*(x- lIF, ) J Fl AG f* i x -VFiY 


and analogously, 


7*(G) — 7&( F j) < / 

JF 


dx 


FjAG f*( x ~ Vg) 
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Since B r n K 1 / 2 = 0, | yp j \ < 1/4, and \yc\ < 1/4, we know that l//*(x — Vf,) > 4/m$ and l//*(x — j/g) > 
4/m$ for any a; £ FjAG , implying that 


- 74-(G')| < 


4 


|F,AG|. 


Therefore, (4.21) becomes 


< $(G) + A 0 \FjAG\, 

where Ao = 8 + 2A + 4/m$, and so Fj is a (Ao, ro)-minimizer for j large enough. 

We now exploit some regularity theorems for sets Fj that are (A, r , 0 )-minimizers 
C 2 set. First, let us introduce a bit of notation. For x £ R", r > 0, and v £ S n ~ l , 


(4.22) 

that converge in L 1 to a 
we define 


C „{x,r) = {y £ R n : \p v (y - x)| < r, | q v {y - x) < r}, 
D v (x,r) = {y £ R n : \p v {y - x)\ < r, \q v {y - ar)| = 0}, 


where q u (y) = y • v and p„(y) = y — (y ■ v)y- We then define the cylindrical excess 
at scale r to be 


exc (E, x, r, v) 


\v E - v\ 


r n— 1 


dn 


n— 1 




of E at x in direction v 


The following regularity theorem for almost minimizers of an elliptic integrand is the translation in the 
language of sets of finite perimeter of a classical result in the theory of currents, see mmmm- For a closer 
statement to ours, see Lemma 3.1 in |13j . 


Theorem 4.5. Let f be a X-elliptic surface tension with corresponding surface energy $. Suppose E is a 
(A, r 0 )-minimizer o/$. For all a <1 there exist constants e and C\ depending on n, A and a such that if 

exc (E, x, r, is) + Ar < e 

then there exists u £ C 1,Q (£>, y (a:, 7')) with u(x) = 0 such that 

C v {x, r/2) n d*E = (Id + ui/)(D ll (x, r/2)), 

IMIc°(D„(a:o 1 r/ 2 )) < Cpr exc(F, x, r, v ) 1,{ - 2n ~ 2) , 
l!V-u|| c o (Diy(a . 0ir/2)) < Ciexc(F,£,7yi/) 1/(2rl_2) , 
and r"[Va] C o, 0(Di/ ( J)r / 2 )) < Ci exc(F, x, r, v) 1/2 . 


Applying Theorem 4.5 as in |8], we come to prove the following statement. 

Theorem 4.6. Let f be \-elliptic with corresponding surface energy 4> and let {Ej} be a sequence of (A, r 0 )~ 
minimizers such that Ej —► E in L 1 , with dE £ C 2 . Then there exist functions tpj £ C 1 (dE ) such that 


dEj = (Id + ifjV E )(dE), 


and WipjWc'tdE) 0. 


Theorem 4.6 implies that we may express dFj as 

dFj = {x + ifjVK ■ x £ dK}, 


where Ht/jllchdA) —> 0. Moreover, barF} = bar A' and \Fj\ = \K\, so Proposition 1.9 and (4.16) imply that 

| 2 


C6$(Fj) > \\ r ‘Pj\\H 1 (dK) — c P<S>(Fj) ■ 


(4.23) 


On the other hand, Fj minimizes Qj , so choosing Ej as a comparison set and using (4.18) and (4.19), we 
have 

4>(F,) < $(£/■) < $(if) + c 2 e 2 < $(AT) + 2c 2 /3$(F,-) 2 . 


By (4.19), P^(Fj) > 0,. Then, using (4.23) and choosing c 2 sufficiently small, we reach a contradiction. □ 
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FIGURE 1. Notation used for I\ and a parallel set E. 


5. The case of crystalline surface tension in dimension 2 

In this section, we prove Theorem |1.5| As in the previous section, we begin by showing the result in a 
special case, and then use a selection principle argument paired with specific regularity properties to reduce 
to this case. 

Let n = 2 and suppose that / is a crystalline surface tension as defined in Definition |1.4[ with $ the 
corresponding anisotropic surface energy. The corresponding Wulff shape I\ C R 2 is a convex polygon with 
normal vectors {vi}f =x . Let us fix some notation to describe K , illustrated in Figure 1. Denote by s, the side 
of K with normal vector choosing the indices such that ,s,; is adjacent to s,; +1 and Sj_i- Let 9 x £ (0, tt) 
be the angle between .s,; and s.j+i, adopting the convention that s n +i = Si. Let Hi be the distance from the 
origin to the side Sj. By construction, 

f{yi) = Hi. (5.1) 

We say that a set E C R 2 is parallel to K if E is an open convex polygon with {v E } = {u, ; }^ 1 , that 
is, ve(x) £ {vi\^L x for all x £ d*E, and for each i £ {l,...,iV}, there exists x £ d*E with v E {x) = 

For a set E that is parallel to K , we denote by cr, the side of E with normal vector v l: and h t the distance 
between the origin and cr,; again see Figure 1. We define e, = hi — Hi. Notice that e* has a sign, with e. ( > 0 
when dist(0, Sj) < dist(0, cu) and a < 0 when dist(0, Si) > dist(0,(7i). For simplicity of notation, we let 
|s| = - H 1 (s) for any line segment s. 

The following proposition proves strong form stability for sets E that are parallel to K such that \E\ = \K\ 
and \EAK\ = inf{\EA(K + y)\ : y £ R 2 }. Then, by a selection principle-type argument and a rigidity result, 
we will reduce to this case. 


Proposition 5.1. Let E C R 2 be parallel to K such that \E\ = \K\ and \EAK\ = inf{|lf A(K+y)\ : y £ R 2 }. 
Then there exists a constant C depending on f such that 

fe{E) 2 < C8*{E). 


Proof. Let E be as in the hypothesis of the proposition. By (|5.l|), we have 

N 

' HArrA. <S>(K\ = V HAsA. IRI=V-^ 


N N h II JV rr I I 

*{E) = Y,Hi kii, m) = Y, H ^ i, i^i = w = E — 


Recalling that Ci = hi — Hi, we may express the volume constraint \E\ = \K\ as 


2 
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FIGURE 2. The surface energy of rK relative to K \E is bounded by the right hand side of (5.51. 


E ^ = m = IU = Y. ^ 


N 

E 


£i\cri\ 

2 


Furthermore, 


N 


N 


2\K\5*(E) = *(E) - $(K) = - |*|) = 


Note that |e*| < C\EAK\ for some constant C = C{f), and so by (1.4), 

/ N \ 2 

<C6*{E), 

' i— 1 ' 

and in particular, |ei| 2 < C5<s>(E) for each i. 

Step 1: We use (2.8) and add and subtract to obtain 


(5.2) 


(5.3) 


2|A| - 2|A| 


M E?<'ME)-f-A-' ) = S , { E) + -L.Xf £-\ 

2|a|V Jej*\ x )' 2|iv| V J k \e f*\ x ) Je\kJ*\ x )' 

Thus we need only to control the term A — B linearly by the deficit, where 

dx „ f dx 


A = 


B = 


Jk\e Je\k f*( x ) 

To bound the term A — B from above, we bound A from above and bound B from below. Our main tool is 


the anisotropic coarea formula in the form given in (2.13). 

wl 

dx 


First, we consider the term A, where (2.13) yields 

$(rK; K \ E) 


A = 


/ K\E 


f*( x ) 


dr= r*(rK;K\E) ir 

Jo r 


(5.4) 


We introduce the notation 


J- = {* e {1, . ■ • N} : c* < 0}, 1+ = {!,.. .N}\I~. 


From (5.4), we obtain an upper bound on A by integrating over r, for each i £ I , the part of the perimeter 


of rK that lies between cu and s,;. This means that for each r, we pick up the part of d*{rK) that is parallel 
to <Ji and Si, as well as part of the adjacent sides: 


§(rK-K\E) < 


u i | , H ( rHi-hi ) o rHi-hi ) 

tlir\Si\ + -H7-1 —;—777-7 h -n ;+\ 


sin(6 l i _i 


sin ( 07 ) 


(5.5) 


see Figure 2 and recall (5.1). This and (5.4) imply that 
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^E 


hi/Hi 


rr , | . rr ( rH i ~ hi) {rH t - hi) 

tlir\Si\ + Mi- 1———---h Mi + \ 


sin(0j_ 1 ) 


sin (0i) 


dr 

7 

r 


(5.6) 


Now we add and subtract the term f h ./ H . The idea is that Hi\(ii\ gives a rough estimate of the 

term in brackets on the right hand side of ( |5.6| ). Indeed, for each r, the part of d* ( rK ) between cr, and Sj has 
length roughly equal to H ,We will see that this estimate is not too rough; the error can be controlled 
by the deficit. Thus we rewrite (|5.6[) as 


E / 

■r- r- J hi [Hi 
l£l ' 


HAaA 


dr + E 


HAsA + 


iei- 


’hi/Hi 


Hi - h -± 
r 




H 


i +1 


sin(0i_i) sin(0j) 


HAaA 


dr. 


Noting that H,/ sin(6j) < C = C(f) for each i,j , the right hand side is bounded by A\ + A 2 , where 


*=E 


iei- 


r 

to 

II 

E 

E 

+ 

0 

\h, - til 

Jhi/Hi r 

V-r_ Jhi/Hi 

l£l 11 1 

r J 


HAoi 


dr. 


dr = 


<S>{rK ; E \ K ) 


The term A 2 is the error term that we will show is controlled by the deficit in Step 2. 

First, we perform an analogous computation for B , and show how, once the error terms are taken care 
of, the proof is complete. Again, by (2.131, we have 

f dx r°° $(rK-,E\K) 

J e\k f*(x) Jo r j ! r 

To bound B from below, we integrate, for each i £ I + , only the part of d*(rK ) that is parallel to S; and cr,; 
and lies between s t and a t . We call this segment t\ := E \ K n {e^ + rxi}, where a is the vector parallel 
to the sides <Xj and Sj, Xi £ Si, and r £ [ 1, hi/Hi}. 


dr. 


Thus, letting s\ be the side of rK parallel to s, ; and recalling ( |5.1[ ), we have 

<&(rK-E\K) 


dr > E 

i£l+ ' 


hi/Hi H‘\s r j nl[| ^ 


Once again, a rough estimate for H^s/ fl l\\ is given by H, ay |. We will again show that this estimate is not 
too rough, specifically, that the error between these integrals is controlled by the deficit. So we continue: 


*>E 

iei +' 


(•hi/Hi 


HAaA 


dr + E 

iei +' 


rhi/Hi HAslnK I HAaA 


dr = Bi + B 2 , 


where 


* = E 

i£l+ ' 


rhi/H * ffikd 


dr, 


B * = E 

ie/+' 


rh< / JTi #*14 n£[| JTiltril 


dr. 


Like A- 2 , B 2 is an error term that we will show is controlled by the deficit in Step 2. 

Before bounding |^4.21 and |A? 2 1 by the deficit, let us see how this will conclude the proof. As we saw, 
/3®(E) 2 < 6$(E) + 2 | 7 f|(^ — #)■ Recalling that hi = Hi + a, 

Bl HAaA , [ hBHi HAaA 


A~B=y] 

iei- Jh i/H, 


dr — 'Y 

rp / ^ 

iei +' 
hi 


-dr + A 2 - B 2 


E lo s (ff) ~ E lo s (^) + a 2 - b 2 

iei~ 1 iel+ 

N N N 

Y. HiWi \+ 0(e ?)) + A 2 - B 2 = - 0(e?) + A 2 - B 2 - 


i=l 


i= 1 


i=l 
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The first term is precisely equal to 2\K\8$(E) by (5.2), while JTO(e?) < C5$(E) by (5.3). Therefore, once 
we show that \A 2 \ and \B 2 \ are controlled linearly by the deficit, our proof is complete. 

Step 2: In this step we bound the error terms. We show that \A 2 \ < C5$(E); the proof that \B 2 \ < C8$(E) 
is analogous. The main idea for estimating the integral A 2 is to show that the contribution of the adjacent 
sides is small, and then estimate the rest of integrand slice by slice. Recalling A 2 , the triangle inequality 
gives 

r 1 w r 1 r u . -i 

(5.7) 


\a 2 \ < 

E [ —(r\si\ - \a. t \)dr 

+^E 

f 

£ 

1 

i | 3- 

dr 



iei- 

Jhi/Hi 

r 



The second term in corresponds to the contribution of adjacent sides. By hi = Hi 


C E 

iei~ 

=C E |- 


hi/Hi 


H r - 


hi 


dr 


CJ2 (Hi - hi) + hi log(E) 


iei~ 


iei~ 


Hr 


+ 0(e 2 i) =CY, + 0(e?) = c£o(e?) < Ctf*(£7). 


iei- 


To bound the first term in (5.7), we will show that |r|si| — |eq|| < <7max{|ej_i|} for r £ [ hi/Hi , 1], where the 
constant C depends on d>, and then obtain our bound by integrating. To this end, we rotate our coordinates 
such that Vi = e 2 , so the side s,; has endpoints (a, Hi) and (b,Hi) for some a < b. We compute explicitly 
the endpoints of cr.;; it has, respectively, left and right endpoints 


(a + tan (6»i_i - tt/2) a - . 2 1 . , hi) 

V sm(6» i _ 1 ) J 


and 


b — tan ( Oi — 7t/2) a + 


e i+i h 1 

siller 


Thus 


o'i = 


b - tan ( 6i - tt/2) a + ~ ( a + tan (cr.^i - tt/2) a - . € ‘ 1 . ) 

sin (Oi) V sm(0j_ 1 )/ 


sin (0*) 

and so 

lki| — \b — a|| < (7(1 ei | + |ej+i| + |ej_i| 
where C depends on $. Therefore, recalling that \b — a\ = |sj|, 


r\Si\ - \<Ji\ 


< (1 - r)\s.i\ + C(]e»] + |e i+ i| + |e*_i|) < ~^~\si\ + C max{| Cj |} < Cmax{|e,j}. 

tin 


Given this estimate on slices, we integrate over r: 


E 

iei~ 


Hi 


'hi/Hi 


rsj — \g. 


i\)dr< C max{|e,-|} V / —dr = Cmax{\ej\} V H t log (-E 
iei~ Jhi / Hi r iei- 


= <7max{|ej|} ^2 ( e i + °( e i)) = 0(max{|e,-| 2 }) < C(^)8^{E), 


ie/~ 


where the last inequality follows from (5.3). 


□ 


We prove Theorem 1.5 after introducing the following definition that we will need in the proof. 
Definition 5.2. A set E is a volume constrained ( e,r]o)-minimizer o/$ if 

$(£) < $(F) + e\EAF\ 

for all F such that |Ej = |F| and (1 — r/o)E C F C (1 + rjo)E. 
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Proof of Theorem m By Proposition |2.4[ we need only to show that there exists some C depending on / 
such that 


M E f < C6§(E). 


(5.8) 


for all sets E of finite perimeter with 0 < \E\ < oo. Suppose for contradiction that (5.8) does not hold. 
There exists a sequence {Fj} such that \Ej\ = \K\, 5<s>(Ej) —» 0, and 


HM < $(F) + c 3 MEj ) 2 


(5.9) 


for C 3 to be chosen at the end of this proof. By an argument identical to the one given in the proof of 
Theorem 1.3 we obtain a new sequence {Fj} with F ;) C B Eo for all j such that the following properties 
hold: 


• each Fj is a minimizer of Qj(E ) = <I>(F) + |/3$(F) 2 —e 2 |+A||F| — \K\ | among all sets E C B Eo , 

where e 3 = /3®(Ej)-, 

• Fj converges in L 1 to a translation of I\\ 

• \Fj\ = \K\ for j sufficiently large; 

• the following lower bound holds for /3$(Fj) : 


e 2 < 2 MEj) 2 - (5.10) 

Translate each Fj such that \FjAK\ = inf{|T)A(A' + y)\ : y £ R 2 }. We claim that for all e > 0, there 
exists ??o > 0 such that Fj is a volume constrained (e, ?yo)-minimizer of $ (Definition |5.2[ ) for j large enough. 
Indeed, fix e > 0 and let rji = Cie, where ci = Ci (/) will be chosen later. By Lemma [2.3[ there exists 772 
such that if (1 — rj 2 )K C E C (1 + 772 )K, then \y E \ < 771 - Let 770 = min{? 7 i, 77 2 }/ 2 . 

By Lemma |3.3[ each Fj satisfies uniform density estimates, and so Lemma |3.4| implies that, for j large, 
(1 — 770 )K C Fj C (1 + rjo)K and thus \yF } \ < Vi- Let E be such that \E\ = |Fj| and (1 — r) 0 )Fj C E C 
(1 + r) 0 )Fj. Then \y E \ < rji and 

(1 - rji)K C Fj C (1 + 7 n )K, (1 - m )K CE C (1 + m )K. 


Because Fj minimizes Qj , 


*( p j) + \M F of - S 2 | < *(E) + ME) 2 ~ e 2 \ 

and so by the triangle inequality and since 777 $ < M$, 

< HE) + ^\ME ) 2 - Mf 3 ) 2 \- 

If <j>(A ? ) < <h( E) 1 then the volume constrained minimality condition holds trivially. Otherwise, with a bound 


as in (3.11), we have 


HFj) < HE) 


*(Fj) - *(E) , \ME)-MH)\ 


2 2 

and so 

$(F j )<$(A) + | 7 $ (A)- 7 $ (F j )|. 

As in the proof of Theorem 1.3 the Holder modulus of continuity for 7 $ shown in Proposition |2.1[ 2) 
does not provide a sharp enough bound on the term 17$(F) — 7 $(Fj)|; we must show that 7 $ is Lipschitz 
when the centers of E and Fj are bounded away from their symmetric difference. In this case, we must be 
more careful and show that the Lipschitz constant is small when |F| = |F| and E and Fj are L°° close. If 
7 $ (F) > 7 $(Fj), then using (2.9), we have 


7$(F)-7$(Fj) < [ 
J E 


dx 


e f*(x- y E ) 


dx 

If, IM^Ve) 


dx 


dx 


' E\Fj f*i x ~ Ve) JFj\E f*i x ~ Ve) 
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One easily shows from the definition that for any x, y £ R ra , 

f*(x) — —\y\ < f*{x -y)< f*{x) + —\y\- 

Therefore, since (1 — r]i)K C EAFj C (1 + rji)K and \y E \ < rji, 

1 - 771(1 + 1/771$) < /*(x -y E )< 1 + ?7i(l + l/m$) 
for a; £ EAFj, implying that 

da; /' da; 


7<&(F) < [ 

Je' 


< Crn\EAFj 


IE\Fj 1 — 7 ?l(l + l/r7T,^>) JFjXE 1 + 7 ? 1 (1 + l/ rn 4>) 

where C = 1 + 1 / 771 $. The analogous argument holds if 7 $ (2?) <7 <s>{Fj), and so 

$(F,) < $(£) + Crft\EAFj\. 

Letting ci == 1/C, we conclude that 7/ is a volume constrained (e, ? 7 o)-minimizer of surface energy, and for j 
large enough, (1 — y 0 /2)K C Fj C (1 + rj 0 /2)K by Lemma |3.4| Therefore, Theorem 5.3 below implies that, 


for j sufficiently large, Fj is a convex polygon with v Ej (x) £ {vi}fL t for "H 1 -a.e. x £ dFj. Moreover, for any 
77, (1 — rf)K C Fj C (1 + r])K for j large enough, so actually {+f. } = {vj}^ for j sufficiently large. In other 


words, for j large enough, Fj is parallel to K , so Proposition 5.1 implies that 


P<s>(Ej) 2 < C\5q>(Fj), 


(5.11) 


where C\ depends on /. On the other hand, Fj minimizes Qj, so comparing against Ej and using (5.9) and 


(5.10) implies 


$(*i) < + c 3 e 2 < &(K) + 2cMFjf. 


□ 


By (5.10), (3${Fj) > 0, so choosing c 3 small enough such that c 3 < \K\/C\, we reach a contradiction. 

Theorem 5.3. TJ1 Theorem 7] Let n = 2 and let f be a crystalline surface tension. There exists a constant 
Co such that if, for some 77 > 0 and some 0 < e < e 3 , (1 — ri/2)K C E C (1 + rj/2)K and E is a volume 
constrained ( e,rj)-minimizer, then E is a convex polygon with 

v E {x) £ forTlf-a.e. x £ dE. 

Remark 5.4. In (T7, Theorem 7], Figalli and Maggi assume that E is a volume constrained (e, 3)-minimizer 
(and actually, their notion of (e,3)-minimality is slightly stronger than ours). However, by adding the 
additional assumption that (1 — rj/2)K C E C (1 + r)/2)K, it suffices to take E to be a volume constrained- 
(e, 77 ) minimizer (with the definition given here) with 77 as small as needed. Indeed, if (1 — r]/2)K C E C 
(l + ? 7 / 2 )A', then (1 — 77 ) 2 ? C co(E) C (1 + 77 ) 2 ? where co(2?) is the convex hull of E. Then, in the proof of [171 
Theorem 7], the only sets F used as comparison sets are such that |2?| = |F| and (1 — if)E C F C (1 + if)E. 

6. Another form of the oscillation index 

The oscillation index /?$(2?) is the natural way to quantify the oscillation of the boundary of a set 2? 
relative to the Wulff shape K for a given surface energy 4>, as it admits the stability inequality (1.8) with 


a power that is independent of /. One may wonder if it would be suitable to quantify the oscillation of 
2? by looking at the Euclidean distance between normal vectors of 2? and corresponding normal vectors of 
K. While such a quantity may be useful in some settings, in this section we show that it does not admit 


a stability result with a power independent of /. This section examines the term /3J(2?) defined in (1.10) 
and gives two examples showing a failure of stability. We then give a relation between /3$ and /3J for 7 -A 
convex surface tensions. As a consequence of Theorem 0 this implies a stability result for /?J, though, as 
the examples show, there is a necessary dependence on the 7 -A convexity of /. 

The following example illustrates that there does not exist a power a such that 

Pi{E) a < C(n, f)6$(E) (6.1) 
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FIGURE 3. The sets Eg are formed by cutting away a zigzag triangle from the top and bottom of 
I\g have 5g{Eg) —> 0 but /3g{Eg) —> oo as 8 —> 0. 


for all sets E of finite perimeter with 0 < \E\ < oo and for all surface energies $. 

Example 6.1. In dimension n = 2, we construct a sequence of Wulff shapes Kg (equivalently, a sequence 
of surface energies $g) and a sequence of sets Eg such that Sg(Eg) -a 0 but fig {Eg) —► oo as 9 — > 0. We use 
the notation Sg = c)$ e and fig = /3§ . 

We let Kg be a unit area rhombus where one pair of opposing vertices has angle 9 < j and the other 
has angle f —9. The length of each side of Kg is proportional to 0 -1 / 2 . Let L = 0 -1//4 . We then construct 
the sets Eg by cutting away a triangle with a zigzag base and with height L from both corners of Kg with 
vertex of angle 9 (see Figure 3). We choose the zigzag so that each edge in the zigzag is parallel to one of 
the adjacent edges of Kg. By taking each segment in the zigzag to be as small as we wish, we may make 
the area of each of the two zigzag triangles arbitrarily close to the area of the triangle with a straight base, 
which is 

A = L 2 tan(0/2) = 6» _1/2 tan(0/2) « 6 » 1/2 , 

as this triangle has base 2Ltan(0/2). Both of the other two sides of the triangle have length m = Lj cos(0/2). 


Let us now compute the deficit Sg and the Euclidean oscillation index /3g of Eg. By construction, $g(Eg) = 
<&g{Kg) = 2, and therefore 

= 2(TY4jIA - 1 = (TT^ - 1 = " ,/2 + 


To compute /3g{Eg) 2 , we cannot characterize the point y for which the minimum in (1.10) is attained in 
general. However, something may be said for an n-symmetric set, i.e., a set E that for which there exist 
n orthogonal hyperplanes such that E is invariant under reflection with respect to each of them. The 
intersection of these orthogonal hyperplanes is called the center of symmetry of E. Indeed, a slight variation 
in the proof of [22 Lemma 5.2] shows that 

3g(E) a Uk-iv-iei 1 /.' 1b 1 - ■ < 6 - 2 > 

where z is the center of symmetry of E. By construction, Eg is a 2-symmetric set with center of symmetry 
0 , so 

9 « (E)2 2 2(1 -A)V* L 1 “ VE - {X) ■ VK -im) m '' 
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FIGURE 4. The sets E r formed by replacing the top and bottom of the i q unit ball with a cone 
show that (6.11 cannot hold for a < 4. 


where Z denotes the union of the two zigzags. By construction, 'H l (Z') is exactly equal to 'H 1 { dKg \ dEg) = 
4m. Moreover, because the edges of Eg are parallel to those of I\g, we find that 


1 - u Eg (x) ■ VK e 


-JL.) = J° 

(1 — COs(7T 


X £ Z 1 

0) x£Z 2 


where Z\ is the set of x £ Z where VE e (x) is equal to vk b ( f ^ x \ ) and Z 2 is the set of a : £ Z where 
v Ee (x) is equal to the normal vector to the other side of Kg. Moreover, we have constructed Eg so that 
U l {Z x ) = n\Z 2 ) = 2to. Thus, as 9 < f, 


j3g{Eg) 2 > ^ / 1 — cos( 7 r — 9) dH. 1 > ^ = m = l/(d 1 ^ 4 cos(0/2)) —» oo 

2 Jz 2 2 

as 9 —> 0. Therefore, for any exponent a, the inequality (6.1| fails to hold; we may choose 9 sufficiently 
small such that Eg is a counterexample. 


The next example shows that even if we restrict our attention to surface energies that are 7 -A convex 
(Definition 1.61, an inequality of the form in (6.1) cannot hold with an exponent smaller than a = 4. The 
example is presented in dimension n = 2 for convenience, though the analogous example in higher dimension 
also holds. 


Example 6.2. Fix p > 2 and define the surface tension f p (x) = (|a?i| p + \x 2 \ p ) 1 ^ p to be the £ p norm in 
R 2 . We show below that f p is a 7- A convex surface tension. Holder’s inequality ensures that the support 
function /* is given by f q , in the notation above, where q is the Holder conjugate of p. The Wulff shape 
K = {f q {x) < 1} is therefore the £ q unit ball. We let denote the surface energy corresponding to the 
surface tension f p . 

We build a sequence of sets {E r } depending on p such that, for any a < 4, we may choose p large enough 
so that S p (E r )//3*(E r ) a —> 0 as r — > 0. Here we use the notation /3*(E ) = /3$ p (E) and S P (E) = 5<s> p (E). 

We may locally parameterize K near (0,1) as the subgraph of the function v q {x\) = (1 — . Thus 

v'qixi) = —\x\\ q ~ 2 Xi/(\ - \xi\ q ) 1/p and 


vk{{x 1 , Vg(Xl))) 


( \x 1 \ q ~ 2 x 1 .\ 

\ (! —I 3 : l| q, ) 1/p ’ j 

yjl + 


(1-\x 1 \i) 2 /p 


{\xi\ q 2 x 1 + Q(\x\\ 2q 1 ), 1 ) 

y/l + |Xi| 2 9 -2 + 0(\Xi\ 3q ~ 2 ) 


(6.3) 


The sets E r are formed by replacing the top and bottom of K with cones. More precisely, let C r = (—r, r) xR. 
We form E r by replacing DK IT C r with the graphs of the functions w and —w, where Wi : (—r, r) —>• R 
is defined by w(xi) = —r 9_1 |a;i|/(l — r 9 ) 1 ^ + Co, with Co = (1 — r q ) x ! q + r q /{ 1 — r q ) 1 / p _ The constant 
Cq is chosen so that w{r) = v q (r) and w(—r ) = v q (—r). For X\ £ (—r,r) for r < 1, we have w'(xi) = 
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—r 9 _ 1 sgn(xi)/(l - r 9 ) 1 ^ and 

v E ((xi, w>(xi))) = 

V ^ (1 -ri) 2 /p 

Now, $ p (E r ) = ® P {K) + $ p (E r ; C r ) - $ P (K; C r ), so 

f r ( r q \Vp / 
$ p (E r )-$ p (K) = J^— + 1 ) -( 


(sgn(xi) ( 1 _ T ,,) 1 / j,, (sgn(xi)r 9 1 + 0(r 2q 1 ), 1 ) 


\/l + r 2q ~ 2 + 0(r 3q ~ 2 ) 


(6.4) 


Pi 


+ 1 


i/p 


dx i 


1 - |Xi|9 

= 1 / r q - \ Xl \ q + 0(r 2q ) dx i = Cr q+1 + o(r 9+1 ). 

P J-r 

The graph of w lies above the graph of v q for all \xi\ < r, so \E r \ > \K\. This implies that 

S r (E r ) < = &«+' + o(r«+ 1 ). 

Next we compute /3*(E r ) in several steps. As in Example 1 6 . 1| E r is a 2-symmetric set with center of 
symmetry 0, thus it is enough to compute the right hand side of ( |6.2[ ). First, the Taylor expansions in (6.31 
and (6.4) imply that, for igC r D d*E , ve{x) • vr () is given by 


(|xi | 9 2 Xi + 0[\x\\ 2q : ), 1 ) (sgn(xi)r 9 1 + 0(r 2q x ), 1 ) 


l + |xi|V ?- 1 + 0 (r 3? - 2 ) 


\Jl + \xi\ 2q ~ 2 + 0 (|xi| 3 i -2 ) yj\ + r 2q ~ 2 + 0(r 3q ~ 2 '> a /(1 + \xi\ 2q ~ 2 + r 2q ~ 2 + 0{r iq ~ i )) 

= 1 + |xi|«-V " 1 - ^(|xi | 29-2 + r 2q ~ 2 ) + O 
For x G d*E \ C r , v E {x) ■ = 0- Hence, 


= 1 + Ix^-V 9 - 1 - ^(\x 1 \ 2q ~ 2 +r 2q ~ 2 ) + 0(r 3q ~ 2 ) = 1 - ^(|xi | 9_1 - r 9 " 1 ) 2 + 0(r 3q ~ 2 ). 


a-E 


: (/*(*)) dn ') 


1/2 


I a-E nc, 


1 — v E ■ vr( ) dH 1 ) 
f* W / 


1/2 


v E (x) - Vr 

= (^J — r q ^ 1 ) 2 \/Y^r 2q ~ 2 ~^~0(r 3q ~ r Ej + 0(r 3q ~ 2 ) dx^j 

= (^j Jd^il 9 " 1 -r 9 ” 1 ) 2 +0{r 3q ~ 2 )dx 1 ' S j 7 =C , P- 1 / 2 + o(r 9 - 1/2 ). 


Furthermore, |A| = |AT| + o(l), so V2\K\ 1 / 4 |A| P 4 = \/2|A'| P 2 + o(l), and so 


P*p(Er) = 


(\L ^ - •*( sfe) l 2 ® 1 ) 1/2 = Cr " ,/2+o(r, ~ ,/2) - 

6 p {E r ) 


v / 2|A'| 1 / 4 |£| 1 / 4 


Therefore, 


P;{E r y 


7*(9+l) 

r a(q-l/2) 


— ? ,9+l -<rq+er/2 


This quantity goes to 0 as r goes to zero if and only if q + 1 — crq + er/2 > 0, or, equivalently, if and only if 
2 (o-\) > d- For any a < 4 we may find 1 < q < 2 (^- 1 ) • Therefore, for any a < 4, there exists a 7 -A convex 
surface tension / such that a bound of the form <5$(A) > C/3^(E) a fails. 


When / is 7 -A convex (recall Definition 1.6), we can control fi%{E) by /3<j>(A). As one expects after the 
previous example, the exponent in this bound depends on the 7 -A convexity of <h. Indeed, this is the content 
of Theorem 1.7 First, we show that the l v norms f p as defined in the previous example are 7 -A convex for 
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each p £ (1, oo). In the case where 1 < p < 2, f p is actually uniformly convex in tangential directions, so it 
is 7-A convex with 7 = 0 . Indeed, f p {y + r) = f p (v ) + V/ p (i/)r + | f* V 2 / p (p + sr)[r, r]ds , and thus 

1 f 1 

f p (v + t) + f p {v - r) - 2f p (v) = - J V 2 f p {u + sr)[r,r]ds. 

We can bound the integrand from below pointwise. We compute 


dufpiy) = ip- 1 )( 


,.\P-2 


\2p-2 


/pM p_1 f P {v) 2p 


-1 / ’ 


-Wph = (1 -p) 1 


IP —2 


"/ a p s 


f P (v) 2p 


-1 


Therefore, if / p (z/) = 1, then 


V 2 / P W = (P-1)E 


2 e i (g)e 1 


2=1 


- (p -!) Z! 

ijJ=l 


h p hih| p 2 p 7 ei(g>e 7 


and so 

n n 2 

V 2 / P (^)[r,r] = (p- 1)^] |Pi| p - 2 r 2 - (p - 1) ( ^ • 

2=1 2—1 

It is enough to consider r such that r is tangent to K p = {/ p < 1} at p, as f p is positive 1-homogeneous 
and the span of v and T„K p is all of K". Observe that V/ p (p) = 5ZT=i \ v i\ p ~ 2v i e i ; this is verified by the 
fact that the support function of f p is f q , and that V/ p (p) = jj^) suc ^ ^at ■ v = f p {v) = \. Thus r 
is tangent to K p at v if and only if r • V/ P (i/) = Xa=i \ v i\ p ~ 2v i T i = 0. Therefore, for such r, 

n 

v 2 / P (^)[r,r] = (p- hl p_27 f ^ (*— 1 )M 2 - 
2=1 

In the case where p > 2, we use Clarkson’s inequality, which states that for p > 2, 

+ , * fx-y\ p ^ fp{x) p , fp{y) p 

J p { 2 J +Jp \ 2 J S 2 + 2 

For v such that / p (p) = 1 and r tangent to K p at v with f p (r) = 1, Clarkson’s inequality with x = v + er 
and y = v — er implies 

2 e p < f p (v + er) p + f p (y - er) - 2 . 

This is almost the condition we need, except we have f p instead of f p for the terms on the right hand 
side. Note that both f p (v + er) and f p (i/ — er) are greater than 1, as moving in the tangent direction to 
K p = {fp < 1} increases f p . The function z p is convex with derivative pz p_1 , so z p < 2 P ~ 1 pz + (2 P ~ 1 p — 1) 
for all z £ [1, 2]. Applying this to z\ = f p (v + er) and z 2 = f p (v — er) yields 

2 e p < 2 P ~ x pf p (y + er) p + 2 p ~ 1 pf p (y - er) p - 2 ( 2 *— x p). 

Thus f p is 7 -A convex with 7 = p — 2 and A = 1/(2 P ~ 2 p). 

The following lemma about 7 -A convexity condition will be used in the proof of Theorem |1.7| 


Lemma 6.3. Assume that f is 7 -A convex. Then for all v,t £ R n such that v ^ 0, 


fiy + r) > 


A 


2 2 +^|i 


( t r) 


v 

v\J \v\ 


2+7 


+ /M + V/(iz) 


(6.5) 


Proof. Note that if / is 7 -A convex, then / is convex. To see that (6.5) holds for given v 0 and r 0 , we let 
f(v) = f{ u ) — f{ v 0 ) — V/( vq) ■ (v — v 0 ). At the midpoint the 7 -A convexity condition gives us the 

following: 


K v 0 ) + f{ v 0 + r o) - 2/(^0 + y) > 


A 

TO 

f T ° 

p 0 > 

^0 

hi 

2 

V 2 

hi 2 

hi 


2+7 


Convexity implies that /(pq + tP) > 0, and /(pq) = 0 by definition of /, implying (6.5). 


□ 
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Finally, we prove Theorem |1.7[ 

Proof of Theorem The quantity (3%{E) measures the overall size of the Cauchy-Schwarz deficit on the 
boundary of E, while /3®(E) measures the overall deficit in the Fenchel inequality. Our aim is to obtain 
a pointwise bound of the Cauchy-Schwarz deficit functional by the Fenchel deficit functional, and then 
integrate over the reduced boundary of E. Without loss of generality, we may assume that \E\ = \K\ = 1 
and E has center zero in the sense defined in Section 12.31 

We fix x £ d*E and consider the Fenchel deficit functional G(v) = f{v) — v ■ which possesses the 

properties that G{v) > 0 and G(v) = 0 i f and only if v = cV/*(x) for some c > 0. 

Let w = |v/l(xj| = Vk (j3x})‘ Lemma 


6.3 


with v = w and t = v E — w, implies that 


A 


fiPE) > xwrdO'E -w)~ ( 0 e - w) ■ w)w | 2+7 + f(w) + Vf(w) ■ (y E - w). 


2 2+ 7 


Therefore, since V/(u>) = and f(w) = V f(w) ■ w, 


GOe) > “ w ) ' w ) w 


2+7 


A 

22+7 

A 

22+7 


(1 -(^. W ) 2 )( 2 + 7)/ 2 

((1 - V E ■ Ml)(l + U E ■ I+)) (2+7)/2 . 


We want to show that there exists some ci such that 

G(y E ) > ci(l - v E • w) (2+7)/2 . 


( 6 . 6 ) 


When wv E > — cq for some fixed 0 < cq < 1, then G(u E ) > — co)( 2+7 ^ 2 (l — v E • ■ u ,)( 2 +tO /2 anc j ( 6.61 

holds. On the other hand, when w ■ v E < — Cq for cq small, we expect that jyyy • v E must also be small and 
so G(v E ) is not too small. Indeed, 


< f(w) = 


w = 


■ cos( 0 i) < cos( 6 >i), 


f*(x) f*(x) 

where 0\ is the angle between w and jqqy- Similarly, 

—c 0 >v E -w = cos( 0 2 ), 

where 0i is the angle between w and v E . Noting that 0 < m$/M$ < 1, so cos - 1 (m$/M$) £ ( 0 , 7 t/ 2 ), we 
let 0o = 2cos _1 (to$/M$) + e, where e > 0 is chosen small enough so that 60 < tt. Letting Co = — cos(0o), 
we deduce that 0i < cos _ 1 (w$/M$) and 0 2 > 60 ■ Then 

x Ixl 


v E < 


cos(0 2 — 0i) < cos (cos 1 (m$/M$) + e) < — M$c e , 


jn— 1 


f*(x) f*(x) 

for a constant c e > 0. Since f{v E ) > m$, we have G{v E ) > M$c e , implying ( | 6 . 6 | ) because 1 — v E ■ w < 2. 
Holder’s inequality and (6.6) imply 

[ 1-VE-w dU 11 - 1 < n n ~ 1 {d*Ey^ 2+ ^ ( [ (1 -v E - w) ( 2 + 7 )/ 2 d'H n " 1 ') 2/(2+7) 

Jd*E ' Jd* E ' 

= cf 2/(2+j) P{E ) 7 /( 2 + 7 ) ( f Ci(l -u E - w) {2+l)/2 dH 

' J a* e 

< c - 2 /( 2 + 7 ) P(F;) 7/(2+7) ( J G(i/ s )d'H"" 1 ) 2/( ‘" +7) 

Dividing by n\K | 1 / n |P| 1 / n and taking the square root, we obtain 

DUE) < + 1/WT) ( n|g fi^ |W 


\ 2 /( 2 + 7 ) 


□ 
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